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VARIABILITY  AND  MIXING  OF  THE  SURFACE  LAYER  IN  THE  TYRRHENIAN  SEA: 

MILEX-80,  FINAL  REPORT 


by 

Jon  Moen 


ABSTRACT 


Data  from  the  1980  Tyrrhenian  Sea  mixed-layer  experiment  (MILEX-80)  are 
reported.  Results  from  the  experiment  are  presented  and  discussed.  The 
high  south-north  gradients  in  the  mixed-layer  depth  are  explained  by  an 
Ekman  pumping  mechanism  (upwelling  and  downwelling)  caused  by  a  maximum  in 
the  west  wind  blowing  through  the  Bonifacio  Strait.  It  is  also  argued  that 
this  band  of  wind  maximum  (or  jet)  is  chiefly  responsible  for  the  pair  of 
counter-rotating  eddies  observed  to  be  quasi -permanent  features  of  this 
zone.  Measured  currents  in  the  area  were  found  to  be  up  to  twice  as  great 

as  the  currents  estimated  from  the  density  field.  This  suggests  a  strong 

barotropic  component  to  the  currents  in  the  north  Tyrrhenian  Sea.  Some 
results  from  mixed-layer  models  are  presented  and  compared  with  the  data. 
These  illustrate  the  need  to  generalize  one-dimensional  mixed-layer  models 
to  3-dimensions  to  properly  model  the  mixed  layer  in  areas,  such  as  the 

MILEX-80  zone,  where  there  is  a  high  spatial  variability  in  the  wind  field. 
The  wind-stress  curl  is  estimated  using  two  independent  methods  and  com¬ 

pared  with  historical  data  for  the  same  area.  The  first  method  involves 
estimating  the  upwelling  and  downwelling  required  to  obtain  a  heat  balance 
in  the  region.  The  second  method  results  from  a  relation  between  inertial 
period  temperature  variance  in  the  thermocline  and  inertial  pulsations  in 
the  Ekman  mass  flux  divergence  caused  by  wind  stress  curl.  The  two  methods 
show  good  agreement  with  each  other  and  with  the  historical  data.  The  same 
divergence  pulsations  are  invoked,  although  less  successfully,  to  explain 
the  observed  inertial  peaks  in  the  anticlockwise  rotary  spectra  below  the 
thermocl i ne. 
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INTRODUCTION 


The  surface  mixed  layer  and  the  highly  stratified  layers  beneath  it 
together  form  a  complex  surface  boundary  layer  to  the  deep  ocean. 

Incoming  energy  from  the  atmosphere  is  initially  stored  in  this  'surface 
layer'  in  the  forms  of  potential  energy,  kinetic  energy,  and  heat  energy. 
These  vary  in  time  due  to  the  local  rates  of  conversion  from  one  form  to 
another  and  due  to  their  local  fluxes  across  both  the  upper  boundary  with 
the  atmosphere  and  the  lower  boundary  with  the  deeper  ocean.  They  vary  in 
space  due  to  spatial  variations  in  the  boundary  fluxes  and  due  to  the  dyna¬ 
mic  response  scales  of  both  the  surface  layer  itself  and  of  the  deeper 
ocean. 

Locally,  and  on  short  time  scales,  it  is  important  to  understand  these 

thermal  and  dynamical  processes  of  the  surface  layer  in  order,  for  example, 
to  improve  the  formulation  of  top  boundary  conditions  for  ocean  circulation 
modelling  (usually  achieved  by  modified  one-dimensional  mixed-layer 
models),  or  to  improve  the  interpretation  of  sea-surface  images  sensed 
remotely  by  satellite.  On  longer  time  scales  these  processes  have  an 
important  effect  on  climate  and  long-range  weather  prediction.  In  addi¬ 

tion,  the  surface  layer  and  its  variability  are  important  in  sound  propaga¬ 
tion  loss  models,  sonar  performance  models,  ambient-noise  models,  and  in 
calculating  the  performances  of  hull -mounted  sonars  and  sonobuoys. 

Several  large  international  experiments  to  study  the  surface  layer  of  the 
ocean  have  taken  place  during  the  past  decade.  Amongst  these  were  the  MILE 
1977  experiment  in  the  North  Pacific  (see  <1>  and  <2>  for  example)  and  the 
JASIN  1978  experiment  in  the  North  Atlantic  <3>.  Research  in  the 

Mediterranean  has  concentrated  on  measurements  taken  from  and  around  the 
French  Bouee  Laboratoire  moored  in  the  Gulf  of  Lions  during  the  1970' s. 
From  this  the  most  useful  data  set,  in  terms  of  mixed-layer  experiments, 
was  that  of  the  COFRASOV  II  expedition  in  1976.  This  data  set,  used  for 

model  testing  by  Klein  <4>,  was  quite  short  (only  10  days  during  July)  and 

no  account  was  taken  of  upwelling  effects.  The  main  purpose  of  Klein's 
model  was  to  demonstrate  the  importance  of  meteorological  variability  in 
time,  whereas  in  the  present  report  we  are  more  concerned  with  spatial 
variability. 

The  surface  layer  in  the  Mediterranean  differs  from  that  of  the  open  ocean 
in  two  important  ways:  firstly,  the  meteorological  forcing  in  the 

Mediterranean  has  a  far  greater  variability  in  space  than  in  the  open  ocean 
and,  secondly,  below  the  seasonal  thermocline  the  vertical  stratitifaction 
is  very  low  in  the  Mediterranean  compared  with  that  in  the  open  ocean. 
Evidence  of  this  can  be  seen  in  the  atlas  by  Robinson,  Bauer  &  Schroeder 
<5>  where,  with  the  exceptions  of  the  Alboran  Sea,  waters  off  the  western 
part  of  the  north  African  coast,  and  regions  of  the  Agean  Sea,  the  ther¬ 
mocline  virtually  disappears  in  the  Mediterranean  during  winter. 

The  high  meteorological  variability,  due  to  the  great  variety  of  wind 
systems  found  around  the  Mediterranean  coastline,  causes  a  similarly  high 
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spatial  variability  in  the  mixed-layer  depth.  Such  variability  arises  as  a 
balance  between  the  local  rate  of  deepening  that,  following  Kraus  &  Turner 
<6>,  is  proportional  to  the  cube  of  the  wind  speed,  and  an  upwelling  or 
Ekman  pumping  mechanism  that  is  proportional  to  the  product  of  the  wind 
speed  and  its  large  scale  (order  of  100  km)  variations.  The  balance,  then, 
is  between  local  and  non-local  (i.e.  large-scale)  forcing  and  the  resulting 
variability,  which  is  confined  entirely  to  the  surface  layer,  is  by  far  the 
most  striking  to  be  found  in  the  Mediterranean.  The  Mediterranean  surface 
layer  thus  takes  on  a  special  importance  that  it  might  not  have,  say,  in 
the  open  ocean  where  a  thermocline  and  generally  greater  vertical  structure 
can  be  found. 
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1  DESCRIPTION  OF  EXPERIMENT 
1.1  Area 


A  mixed-layer  experiment  (hereafter  referred  to  as  MILEX-80)  was  conducted 
in  the  north  Tyrrhenian  Sea  in  October  1980.  Its  main  purpose  was  to 
investigate  spatial  and  temporal  variability  of  the  surface  layer  during 
the  period  of  transition  from  summer  to  winter  conditions.  The  area  of  the 
experiment,  shown  in  Fig.  1,  was  an  approximately  200  km  square  zone 
laying  to  the  east  of  Corsica,  the  Strait  of  Bonifacio,  and  Sardinia. 

Previous  measurements  in  this  area  <7>  had  indicated  high  variability  in 
the  mixed-layer  depth,  with  a  relatively  shallow,  cold  mixed  layer  to  the 
north  of  the  zone  and  a  relatively  deep,  warm  mixed  layer  to  the  south. 

Satellite  infrared  images,  such  as  the  one  shown  in  Fig.  2  indicate  the 
presence  during  the  summer  and  autumn  of  a  large  cold  patch  of  water  to  the 
east  of  Corsica  and  the  Strait  of  Bonifacio.  In  winter  the  cold  patch 
remains  but  is  less  well  defined  due  both  to  the  merging  of  its  northern 

and  eastern  boundaries  with  the  extended  zone  of  coastal ly  up-welled  water 
off  the  west  of  Italy  and  to  surface  cooling  of  all  the  Mediterranean 
waters.  The  persistence  of  the  frontal  zone  separating  the  cold  water  from 
the  warmer  water  to  the  south  has  been  studied  using  satellite  infrared 
images  by  Philippe  &  Harang  <8>. 

Dynamic  topography  estimates  by  Krivosheya  &  Ovchinnikov  <9>,  using  data 
collected  by  French  and  Italian  hydrographic  research  vessels  during  the 
1950' s,  indicate  the  presence  of  a  permanent  cyclonic  eddy  to  the  west  and 
north  of  the  Bonifacio  Strait,  corresponding  to  the  region  of  cold  water, 
and  a  smaller,  less  intense  anti  cyclonic  eddy  in  the  warmer  water  to  the 
south.  The  authors  attributed  these  eddies  to  forcing  by  the  wind-stress 

curl  (estimated  from  mean  geostrophic  wind)  arising  from  a  band  of  maximum 

of  the  northwest  wind  through  the  Bonifacio  Strait  during  the  summer  months 
July  to  September.  However,  no  explanation  was  given  as  to  why  these 
eddies  should  remain  throughout  the  winter  in  the  apparent  absence  of  this 
wind  maximum. 

A  recent  analysis  of  wind-stress  patterns  over  the  Mediterranean  has  been 
made  by  May  <10>  in  which  surface  stresses  calculated  from  wind  vector 
observations  by  shipping  were  averaged  into  0.2°  squares  for  each  month  of 
the  year.  These  data,  approximately  one  million  observations  covering  the 
20  years  from  1950  to  1970,  indicate  a  band  of  wind-stress  maximum  directed 
more  eastward  from  the  Strait  of  Bonifacio  than  that  cited  by  Krivosheya  & 
Ovchinnikov  <9>.  The  data  also  show  the  persistence  of  the  band  throughout 
the  year. 


1.2  Method 

The  method  of  the  experiment  was  to  monitor  the  deepening  of  the  mixed 
layer  and  erosion  of  the  thermocline  using  an  arrangement  of  taut  moored 
thermistor  chains  and  current  meters.  Six  mooring  sites  A,  B,  C,  D,  E  &  F 
were  established  at  the  positions  indicated  in  Fig.  3.  These  positions  can 
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FIG 


1  AREA  OF  MILEX  80 
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FIG.  2  SATELLITE  INFRARED  IMAGE  FOR  27  OCT,  1980,  SUPPLIED  BY  THE 
UNIVERSITY  OF  DUNDEE.  Light  shades  represent  colder  temperature. 
Letters  A-F  indicate  the  mooring  positions  and  the  line  the  track 
used  for  vertical  sections . 


FIG.  3  POSITIONS  OF  CTD  STATIONS  RELATIVE  TO  MOORING  POSITIONS. 
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also  be  seen  relative  to  the  cold  water  patch  in  the  satellite  image  of 
Fig.  2.  Each  mooring  consisted  of  a  Bendix  surface  float  anchored  by 
approximately  100  m  of  3  trm  steel  cable  for  attachment  of  instruments,  a 
variable  length  (depending  on  water  depth)  of  stretched  6  mm  nylon  Samson 
rope,  and  approximately  100  m  of  3  mm  steel  cable  attached  in  line  with 
buoyancy  spheres,  an  acoustic  releaser,  and  an  anchor  (see  Fig.  4).  One 
Aanderaa  thermistor  chain  of  11  thermistors  was  attached  to  the  upper  steel 
cable  of  each  mooring  and  adjusted,  on  site,  to  span  the 
mi xed-1 ayer/upper-thermocl ine  region.  The  depth  of  each  thermistor  rela¬ 
tive  to  the  uppermost  was  9.7,  12.7,  15.7,  18.7,  21.7,  24.7,  27.7,  32.7, 
37.7  and  42.7  m.  The  sampling  interval  was  set  to  10  min.  One  NBA  current 
meter  was  attached  in-line  under  each  thermistor  chain  and  set  to  sample  at 
5-min  intervals. 

Two  additional  taut  moorings,  C2  and  C3,  each  of  two  current  meters,  were 
deployed  at  position  C  (see  Fig.  4).  One  of  these,  C3,  consisted  of  a 
Neil-Brown  3-axis  acoustic  current  meter  (sampling  interval  1  min  at  4  m) 
and  an  EG&G  VACM  (vector-averaging)  current  meter,  sampling  interval  1.875 
min  at  14  m.  The  other  was  comprised  of  VACM  current  meters  at  23  m  and 
32  m,  both  sampling  at  1.875  min.  The  total  of  five  current  meters  at 
position  C  were  arranged  to  span  the  mixed  layer  and  upper  thermocline  with 
the  aim  of  measuring  vertical  shear  during  deepening  periods. 

Wind  speed,  wind  direction,  air  temperature,  sea-surface  temperature,  and 
barometric  pressure  were  measured  by  an  Anderaa  meteo-chain  mounted  at  a 
height  of  3.5  m  on  the  surface  float  of  mooring  A  (Fig.  4).  The  recording 
periods  for  the  instruments,  together  with  their  depths,  sampling  inter¬ 
vals,  and  mooring  identification  are  given  in  Fig.  5. 


1.3  Onboard  Measurements 

The  moorings  were  deployed  and  recovered  using  the  SACLANTCEN  research 
vessel  MARIA  PAOLINA  G.  (MPG).  During  the  deployment  phase  19  CTD  stations 
were  taken  to  a  depth  of  500  m  and  162  XBT  casts  made  to  300  m.  During  the 
recovery  phase  17  CTD  stations  and  88  XBT  casts  were  made.  An  additional 
52  XBT  casts  were  made  by  the  Italian  Navy  hydrographic  research  vessel 
MIRTO  during  an  intermediary  phase  that  included  surveillance  of  the 
moorings.  All  CTD  positions  are  shown  in  Fig.  3.  The  positions  of  the  XBT 
measurements,  which  were  taken  along  the  ship's  track  between  mooring  posi¬ 
tions  and  CTD  stations,  are  shown  in  Fig.  6. 

Wind  speed  and  direction  at  10  m  were  recorded  from  the  ship's  meteo-chain 
and  sea-surface  (4  m  depth)  temperature  (SST)  was  recorded  by  a  thermistor 
projecting  down  through  a  tubular  opening  near  the  ship's  well. 

The  measurement  periods  for  the  CTD,  XBT,  and  SST  measurements  are  shown  in 
Fig.  5.  Also  included  in  that  figure,  as  an  indication  of  the  weather  con¬ 
ditions  during  the  experiment,  are  the  wind  speed  and  direction  measured  by 
ship's  anemometer  until  13  October  and  thereafter  by  the  meteo-chain  at 
mooring  A.  It  can  be  seen  that  at  least  three  storms  occurred  during  the 
experiment,  in  which  the  wind  speed  was  greater  than  10  m/s.  These  were 
roughly  during  the  periods  9-11,  16-19,  and  25-27  October. 
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Also  available  from  the  MIRTO  ship's  log,  at  approximately  1  h  intervals, 
are  wind  speed  &  direction,  barometric  pressure,  and  wet  &  dry  air  tem¬ 
perature. 


1.4  Data  Recording  and  Preparation 

Data  from  the  moored  instruments  were  recorded  on  magnetic-tape  cassettes, 
transcribed  onto  8-track  HP  tapes,  transferred  to  UNI  VAC  disc  storage,  and 
then  prepared  for  analysis  by  the  replacement  of  'bad'  data  values  and 
truncation. 

Onboard  acquisition  of  XBT  &  CTD  data  was  by  use  of  SACLANTCEN's  lOS  system 
(integrated  oceanographic  system)  on  a  Hewlett  Packard  21MX  computer.  Data 
were  stored  on  HP  8-track  tapes  and  later  transferred  to  UNIVAC  1106  disc 
storage  where  they  were  cleaned  by  the  removal  of  spikes  and  bad  profiles. 

Apparent  wind  speeds  &  directions  and  sea-surface  (4  m)  temperatures  were 
recorded  on  strip  charts  and  later  digitized  and  corrected  before  storage 
on  disc. 

The  XBT  data  collected  by  the  Italian  Navy  were  also  digitized  and  trans¬ 
ferred  to  UNIVAC  disc  storage. 


aanderaa  meteo  chain 


FIG.  4  SCHEMATIC  REPRESENTATION  OF  MOORING  DESIGN. 
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FIG.  5  DEPTH,  SAMPLING  INTERVAL,  AND  RECORDING  PERIODS  FOR  EACH  OF  THE 

INSTRUMENTS.  Wind  speed  indicates  weather  conditions  during  the 
various  phases  of  MILEX  80. 

XBT  TRACKS  MILEX  80 


FIG.  6  POSITIONS  OF  XBT  CASTS 

(a)  Phase  I,  (b)  Phase  II  and  (c)  Phase  III. 
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2  ANALYSIS  AND  INTERPRETATION:  SPATIAL  DESCRIPTION 

2.1  Vertical  Sections 

The  selected  vertical  sections  of  temperature,  salinity,  density,  and  sound 
speed  shown  contoured  in  Figs.  7,  8  and  9  illustrate  the  spatial  charac¬ 
teristics  of  these  variables  in  the  upper  layer  in  the  MILEX-80  zone. 

Figure  7a  is  a  south-north  vertical  section  of  temperature  calculated  from 
XBT  data  taken  between  12  and  13  Oct.  It  shows  a  strong  south-north  gra¬ 
dient  in  the  mixed  layer — being  relatively  shallow  and  cold  (approx.  15  m 
&  19.5°C)  to  the  north  and  relatively  deep  and  warm  (approx.  50  m  &  21.5°C) 
to  the  south.  The  same  general  structure  remains  (Fig.  7b)  after  the 
passage  of  two  fairly  intense  storms  of  several  days  duration  and  wind 

speeds  in  excess  of  10  m/s.  The  new  structure  shows  that,  in  the  cold 

water  region  to  the  north,  the  mixed  layer  has  deepened  by  about  13  m  and 
cooled  by  about  2.5°C,  whereas  in  the  warm  water  to  the  south  the 

deepening  was  about  10  m  and  the  cooling  only  1°C. 

It  is  a  good  measure  of  the  spatial  variability  that  the  south-north  dif¬ 
ference  in  mixed-layer  depth  is  two  or  three  times  greater  than  the  change 
in  layer  depth  at  either  location  caused  by  two  weeks  of  intensive  wind 

mixing. 

Based  on  the  respective  depths  of  the  two  mixed  layers,  the  difference  in 
deepening  rate  between  north  and  south  is  much  less  than  might  be  expected. 
Shallow  mixed  layers  deepen  more  rapidly  than  deep  ones,  because  for 
shallow  mixed  layers  a  greater  fraction  of  the  turbulent  'mixing'  energy, 
created  near  the  surface,  arrives  at  the  mixed-layer/  thermocl ine-interface 
zone.  Furthermore,  the  wind-driven  inertially  oscillating  bulk  (or  slab) 
flow  of  a  shallow  mixed  layer  is  far  greater  than  that  of  a  deep  one,  and 
this  means  that  a  relatively  larger  production  of  turbulent  mixing  energy 
can  be  expected  in  the  mixed-layer/thermocl ine-shear  zone  for  shallow  mixed 
layers  <11>. 

Better  estimates  for  deepening  and  cooling,  in  which  high-frequency 
variations  in  time  have  not  been  smoothed,  can  be  determined  from  the  ther¬ 
mistor  chain  data  presented  in  the  next  chapter.  These  indicate  a  general 
deepening  and  cooling  of  the  mixed  layer  in  the  north  of  about  11.5  m  and 
2.1°C  respectively  and  in  the  south  of  about  6  m  and  1°C  respectively. 
Although  this  represents  a  greater  relative  increase  in  depth  for  the  more 
shallow  mixed  layer  to  the  north  than  for  the  deeper  one  to  the  south,  the 
north-south  ratio  for  deepening  rates  is  still  far  too  low. 

Let  this  ratio  be  denoted  by  R.  Then,  using  the  above  values  of  11.5  m  and 
6  m,  we  obtain  the  value  R=2  for  the  2-week  period.  If  now  we  consider  the 
deepening  rate  due  to  the  downward  penetration  of  turbulence  created  near 
the  surface  <6>: 


dh  u3 

dt  “  ’  (Eq.  1) 
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where  U  is  the  wind  speed,  h  the  depth  of  the  mixed  layer,  and  r  the 
thermocline  gradient  below  the  jump  zone.  We  can  write: 


h2  U3  r2 
2  1 


(Eq.  2) 


where  the  subscript  i  refers  to  the  north,  2  to  the  south  and  we  have 
assumed  the  wind  speed  at  the  two  locations  to  be  the  same.  Taking 
hi  =  25  m  and  h2  =  50  m,  Eq.  2  gives  a  value  of  R  =  4.  Note  that  we 
have  not  considered  the  mechanism  due  to  Pollard,  Rhines  &  Thompson  <11> 

mentioned  earlier,  which  may  in  fact  be  more  important  than  the  Kraus  & 
Turner  mechanism  <6>  for  the  shallow  mixed  layer  to  the  north.  In  this 
case  the  value  of  R  =  4  would  be  an  underestimate  of  the  deepening  rates 

ratio.  It  is,  however,  sufficiently  high,  compared  with  the  measured  value 
of  R  =  2,  to  indicate  that  some  process  other  than  mixing  is  important  in 

the  MILEX-80  zone.  This  process,  which  will  be  developed  more  fully  in  a 

later  chapter  on  Ekman  pumping,  is  clearly  one  of  upwelling  to  the  north 
and/or  downwelling  to  the  south.  In  view  of  the  presence  of  a  large  cold 
patch  observed  in  the  area  throughout  the  year,  the  emphasis  should  be  on 
upwelling. 


Two  further  examples  of  the  south-north  and  southeast-northwest  structure 
can  be  seen  in  the  vertical  sections  of  Fig.  8.  These  are  taken  from  XBT 
casts  taken  during  phase  II  of  MILEX-80  when  the  MIRTO  steamed  south  and 
southeast  from  A  to  F  and  then  returned  along  approximately  the  same  course 
towards  A.  The  two  sections  are  thus  separated  in  time  by  a  maximum  to  the 
north  of  one  day  and  it  can  be  seen  that  during  this  period  the  cold  water 
has  advected  southward,  apparently  sharpening  the  frontal  boundary 
separating  north  from  south.  The  southern  portion  of  the  section  in 
Fig.  8b,  however,  contains  some  east-west  structure,  with  warmer  water 
lying  to  the  west  of  F.  It  would  appear  that  during  this  period, 
14-15  Oct,  cold  water  was  advecting  in  from  the  north  or  northwest  and  this 
can  explain  why  the  mixed  layer  at  F  has  decreased  from  the  50  to  60  m 
level  seen  on  both  12  and  24  Oct,  shown  in  Fig.  7,  to  a  depth  of  only  35  m. 
This  advection  event  can  also  be  seen  on  the  thermistor  chain  data 
displayed  in  the  next  chapter,  where  it  will  be  discussed  in  more  detail. 

Vertical  south-north  sections  of  salinity,  density,  and  sound  speed  for 
phase  I  are  shown  in  Fig.  9a,b,c.  These  are  less  accurate  than  the  tem¬ 
perature  sections,  both  because  fewer  CTD  stations  were  available  than  XBT 
casts  and  because  the  data  were  spread  over  one  week  in  which  an  appre¬ 
ciable  wind  mixing  and  cooling  took  place.  Nevertheless,  they  give  a  good 
indication  of  the  spatial  gradients  involved.  The  change  in  sound  speed  is 
striking  —  changing  from  1530  m/s  at  the  50  m  level  in  the  south  to 
1512  m/s  in  the  north.  Such  an  18  m/s  change  over  90  km  could  have  a 
significant  effect  on  acoustic  propagation  losses  and  ambient-noise  models. 


2.2  Horizontal  Sections 

Dynamic  anomalies  relative  to  the  450  m  depth  level  were  calculated  for 
each  CTD  station  by  the  usual  method  of  transformation  from  z  to  pressure 
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FIG.  7  VERTICAL  SOUTH-NORTH  TEMPERATURE  SECTIONS  ALONG  TRACK  F,  E,  C,  A. 
(a)  12-13  Oct,  1980,  (h)  23-29  Oct,  1980. 

XBT  cast  mmbers  are  indicated  along  the  x-axis. 
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FIG.  8  VERTICAL  TEMPERATURE  SECTION  FROM  XBT  CASTS  TAKEN  BY  ITALIAN  NAVY 
HYDROGRAPHIC  RESEARCH  VESSEL  ’MIRTO ’ . 

(a)  14-15  Oct,  1980,  (b)  15  Oct,  1980. 

XBT  cast  numbers  are  indicated  along  the  x-axis. 
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VERTICAL  CONTOURS  OF  SOUND 


FIG.  9  VERTICAL  SECTIONS  OF  (a)  SALINITY,  (b)  SIGMA-t  AND  (c)  SOUND 
SPEED  ESTIMATED  FROM  CTD  STATIONS  MADE  BETWEEN  6  AND  13  OCT,  1980. 
CTD  station  numbers  are  indicated  along  the  x~axis. 
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co-ordinates  and  integration  of  specific  volume  anomaly,  6  =  a(T,S,p)- 

a(0,35,p),  between  the  reference  pressure  of  450  dbar  and  the  desired 
pressure  surface.  Contours  of  dynamic  anomaly,  in  dynamic  metres,  for  the 
20  dbar  surface  are  shown  in  Fig.  10  in  which  Fig.  10a  represents  phase  I, 

the  mooring  deployment  phase  of  the  experiment,  and  Fig.  10b  represents 

phase  III,  the  mooring  recovery  phase.  These  show  a  cyclonic  gyre 

underlying  the  cold  water  patch  and  evidence  of  an  anti  cyclonic  motion  to 
the  south.  The  results  are  thus  consistent  with  those  of  Krivosheya  and 

Ovchinnikov  <9>. 

Geostrophic  currents  calculated  from  the  dynamic  anomaly  estimates  were 

found  to  agree  well  in  direction  with  measured  currents  (to  be  illustrated 
later)  but  to  be  about  half  the  speed.  This  is  an  interesting  result, 

since  although  some  discrepancy  between  the  actual  current  and  that  calcu¬ 
lated  from  the  density  field  is  expected  in  the  Mediterranean,  no  measure¬ 
ments  have  been  reported  in  the  literature.  Also,  that  this  discrepancy 
should  be  so  high  is  surprising.  Further  discussion  of  this  result  will 
appear  in  Ch.  4. 

Comparing  Figs.  10a  and  b  we  see  that  the  gyre  appears  to  have  either  moved 
or  to  have  partially  broken  up  during  the  two  weeks  of  stormy  weather  bet¬ 
ween  measurements.  An  insufficient  spatial  coverage  of  CTD  stations  pre¬ 
vents  this  ambiguity  from  being  resolved. 

Figures  11  to  16  present  horizontal  contours  of  the  temperature,  salinity, 
density,  and  sound  speed  at  different  depths.  The  salinity,  density  and 
sound  speed  contours  of  Figs.  13,  14  and  15,  respectively,  were  estimated 
from  the  CTD  data  and  are  for  the  depths  20  m  (or  30  m)  &  100  m.  The  tem¬ 
perature  contours  of  Figs.  11  and  12  are  from  greater  coverage  provided  by 

the  XBT  data  from  phases  I  and  III  respectively;  they  show  that  the  anoma¬ 
lies  in  temperature,  which  form  the  major  contribution  to  the  density 
field,  have  virtually  disappeared  below  100  m.  The  density  contribution  to 
the  geostrophic  flow  is  thus  confined  to  a  relatively  thin  surface  layer 
and  it  is  this  layer  that,  for  most  of  the  Mediterranean,  contains  the  bulk 
of  the  variability  in  the  state  variables  and  in  sound  speed.  Figure  16 
shows  that  the  anomalies  in  salinity  go  somewhat  deeper  than  those  of  tem¬ 
perature  but  that  they  also  effectively  disappear  by  300  m. 


2.3  SST  profiles 

Figure  17  plots  the  sea-surface  (4  m)  temperatures  recorded  along  the  hori¬ 
zontal  tracks  indicated  in  Fig.  18. 

These  records  further  illustrate  the  general  impressions  of  a  south-north 
temperature  change  of  about  3°C,  as  shown  in  profiles  3,  5  &  15,  and  of  a 
cooling  of  this  surface  water  throughout  the  experiment  by  about  3°C  as 
seen  by  comparing  profiles  1  to  6  with  7  to  15.  They  also  illustrate  a 
strong  frontal  zone  separating  the  warm  and  cold  water,  as  seen  in  profiles 
5,  7,  13  and,  especially,  15,  which  shows  a  change  in  temperature  of  2°C  in 
about  2  km.  Such  fronts  are  also  evident  after  quantitative  analysis  of 
satellite  infrared  images  of  the  area  <12>  but  these  will  not  be  presented 
here. 
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FIG.  10  20  m  DYNAMIC  ANOMALIES  CONTOURS. 

(a)  6-14  Oct,  1980,  (b)  23-30  Oct,  1980. 


North 


(a) 

6/14  0ct.80 


Cb) 

North 


20m 


100m 


FIG.  11  TEMPERATURE  CONTOURS  FOR  6-14  Oct,  1980, 
(a)  20  m,  (b)  100  m. 
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FIG.  12  TEMPERATURE  CONTOURS  FOR  23-30  Oct,  1980. 
(a)  20  m,  (b)  100  m. 
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FIG.  13  SALINITY  CONTOURS  FDR  6-14  Oct,  1980. 
(a)  20  m,  (b)  100  m. 
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FIG.  14  DENSITY  (sigma-t)  CONTOURS  FOR  6-14  Oct,  1980. 
(a)  30  m,  (h)  100  m. 


(a) 

(b) 

North 

6/14  Oct,80 

North 

FIG.  15  SOUND  SPEED  CONTOURS  FOR  6-14  Oct,  1980. 
(a)  20  m,  (b)  100  m. 
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FIG.  16  SALINITY  OONTOURS  FOR  6-14  Oct,  1980. 
(a)  200  m,  (b)  300  m. 


19 


SACLANTCEN  SR-75 


•  Day  (Octobar  19aO>-Tima  (QMT) 


FIG. 
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SEA-SURFACE  TEMPERATURE  (SST)  PROFILES. 
Numbers  represent  tracks  indicated  in  Fig. 


18. 


FIG.  18  TRACKS  OF  SEA-SURFACE  TEMPERATURE  PROFILES  PLOTTED  IN  FIG.  17. 
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A  theory  for  the  formation  of  fronts  in  areas  of  non-zero  wind-stress  curl 
and  weak  coupling  of  the  mixed  layer  with  the  atmosphere  can  be  found  in 
Welander  <13>.  The  theory  was  applied  to  large  (~  1000  km)  scales  but  if 
one  can  assume  a  sufficiently  weak  coupling  of  the  mixed  layer  with 
atmosphere  at  smaller  scale  (e.g.  cooling  greater  in  warm  water  region)  the 
mechanism  might  also  be  applied  in  the  MI  LEX-80  zone. 


2.4  Current  Measurements 

Currents  were  measured  approximately  below  the  thermocline  (~  70  m  level) 
at  each  of  the  mooring  positions.  Other  current  measurements  were  also 
made  at  mooring  position  C,  but  since  currents  will  also  be  discussed  as 
temporal  phenomena  in  Ch.  3  these  will  not  be  illustrated  here.  The 
measurements  to  be  presented  here  give  a  spatial  impression  and  have  been 
included  for  that  reason. 

Figures  19  and  20  are  progressive  vector  diagrams  for  each  of  the  mooring 
positions,  with  each  of  the  eight  maps  representing  a  2-day  period  and  each 
progressive  vector  sub-divided  into  12  h  intervals.  Figure  19  covers  the 
period  13-20  Oct  and  Fig.  20  covers  21-28  Oct.  This  representation  thus 
gives  both  spatial  and  temporal  information  for  the  MILEX-80  zone.  The 
directions  of  the  currents  are  consistent  with  the  eddy  structure  indicated 
earlier  by  the  surface  contours  of  dynamic  anomaly  —  a  cyclonic  eddy  to  the 
north  and  an  anti  cyclonic  one  to  the  south.  The  magnitude  for  several  of 
the  2-day  periods  approaches  50  cm/s,  which  is  considerably  higher  than  the 
estimated  geostrophic  currents.  This  observation  was  tested  more  carefully 
by  comparing  the  geostrophic  current  relative  to  450  m,  computed  from  the 
CTO  profiles  surrounding  each  mooring  position,  with  values  of  the  measured 
current,  when  available.  Since  baroclinic  adjustment  can  take  in  the  order 
of  a  day  or  more,  the  measured  currents  were  low-pass  filtered  to  60  h 
before  comparison  was  made  in  order  to  make  them  more  representative  of  the 
density  currents.  Table  1  shows  the  comparison,  together  with  additional 
information  such  as  lag  time  between  CTD  and  current  measurement  (often  CTD 
stations  were  taken  at  the  mooring  position  before  deployment  or  after 
recovery  of  the  current  meter),  type  of  current  meter,  and  depth  of  current 
meter.  There  appears  to  be  reasonably  good  agreement  in  direction, 
generally  within  about  20°  except  for  position  D  which  showed  a  47°  dif¬ 
ference,  but  the  speeds  of  the  measured  currents  are  about  twice  as  high  as 
those  of  the  estimated  currents. 


2.5  Errors  in  current  estimates  and  measurements 

Errors  in  the  estimated  geostrophic  current  may  occ.ur  from  the  choice  of  a 
rather  shallow  reference  level  of  450  m.  Geostrophic  estimates  by 
Krivosheya  &  Ovchinnikov  <9>  for  the  same  area  as  MILEX-80,  but  using  a 
reference  level  of  1000  m,  indicate  that  the  flow  at  450  m  can  be  as  high 
as  25  %  of  that  at  the  surface. 

Another  source  of  error  in  the  comparison  is  in  the  use  of  surface-moored 
NBA  current  meters.  This  error  arises  because  of  the  vertical  oscillations 
of  the  mooring  line  to  which  the  current  meter  is  attached;  these  can  be 
quite  large  when  seas  are  high.  Since  the  drag  on  the  current  meter  is  not 
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FIG.  19  FOUR  2-DAY  PROGRESSIVE  VECTOR  DIAGRAMS  FOR  13-20  Oct,  1980. 
Divisions  on  vectors  are  at  12  h  intervals. 
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FJG.  20  FOUR  2-DAY  PROGRESSIVE  VECTOR  DIAGRAMS  FOR  21-28  Oct ,  1980. 
Divisions  on  vectors  are  at  12  h  intervals . 
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symmetric  (higher  drag  on  the  tail  fins  than  on  the  nose)  and  since  the 
current  meter  is  free  to  rotate  ±  10°  in  the  vertical  plane,  these  vertical 
oscillations  tend  always  to  increase  the  rotor  rotations  and  create  an 
overestimation  of  the  mean  current. 

The  magnitude  of  this  'pumping'  error  is  difficult  to  estimate  since  we 
know  neither  the  coupling  of  current  meter  tilt  with  vertical  oscillations 
nor  the  cosine  response  of  the  NBA  rotor.  Furthermore,  the  amplitude  and 
frequency  of  the  vertical  oscillations  are  not  known,  since  the  surface 
waves  which  they  are  forced  to  follow  (surface  float  and  current  meter  con¬ 
nected  by  steel  cable)  were  not  measured.  We  can  only  note  from  Table  1 
that  the  VACM  current  meters,  which  are  less  prove  to  pumping,  indicate 
currents  of  the  same  order  or  higher  than  the  NBA  current  meters.  The 
current  measured  by  the  VACM  at  32  m  on  8  Oct,  for  example,  was  twice  that 
of  the  NBA  at  54  m.  Both  current  meters  were  below  the  Ekman  layer  during 
this  period  and  so  should  not  be  measuring  the  wind-driven  current.  In  any 
case,  the  winds  during  this  period,  as  during  other  periods  when  CTD 
measurements  were  taken,  were  very  light  and  the  pumping  of  the  NBA  current 
meters  would  therefore  have  been  low.  The  NBA569  and  VACM357  current 
meters  can  be  compared  by  using  the  stick  diagrams  presented  in  Sect.  3.2. 
Since  the  NBA  should  be  measuring  consistently  lower  than  the  VACM  (which 
is  near  the  surface)  but  is  in  fact  measuring  slightly  higher  we  can  deduce 
that  it  is  indeed  overestimating  the  current. 

Other  comparisons  between  measured  currents  and  geostrophic  estimates  will 
be  made  for  the  MILEX-80  zone  and  the  Balearic  sea,  but  these  will  be  pre¬ 
sented  elsewhere. 

TABLE  1 

COMPARISON  OF  ESTIMATED  GEOSTROPHIC  CURRENTS  WITH  MEASURED  CURRENTS 

Lag  time  indicates  the  time  difference  between 
measured  (filtered)  current  and  estimated  geostrophic  current 


MEASURED 

CURRENT 

GEOSTROPHIC  ESTIMATE 

TIME  LAG 
(hrs) 

MOORING 

CURRENT 

METER 

DEPTH 

(m) 

MAGNITUDE 

(cra/s) 

DIRECTION 

(degrees) 

MAGNITUDE 

(cm/s) 

DIRECTION 

(degrees) 

C 

VACM 

14 

33-8 

163.2 

16.5 

163.4 

-12 

VACM 

23 

24.8 

152.6 

13.5 

167.2 

-17 

VACM 

32 

30.5 

178.3 

11.2 

169.8 

-17 

NBA 

54 

14-9 

177.7 

7.5 

173.1 

0 

B 

NBA 

63 

13-1 

324.0 

6.3 

339.6 

0 

D 

NBA 

60 

11.2 

83.8 

3.1 

61.2 

0 

E 

NBA 

68 

6.4 

65.0 

10.2 

18.3 

+  7 

C 

VACM 

14 

24.7 

142.6 

14.2 

132.5 

+  72 

VACM 

23 

14.7 

157.2 

12.8 

124.0 

+36 

VACM 

32 

17.1 

149.4 

10.5 

130.4 

+36 

NBA 

54 

16.0 

152.0 

5.7 

125.1 

+36 

24 


SACLANTCEN  SR-75 


3  ANALYSIS  AND  INTERPRETATION:  TEMPORAL  DESCRIPTION 
3.1  Thermistor  chain  data 

In  response  to  a  number  of  passing  storms  that  occurred  during  the  three 
weeks  of  MILEX-80,  the  mixed  layer  deepened  by  about  15  m  and  cooled  by 
about  3°C.  An  indication  of  the  variability  during  this  period  is  shown  by 
Fig.  21,  which  plots  the  depth  of  the  maximum  thermocline  gradient  as  esti¬ 
mated  from  the  thermistor  chain  data.  Effectively  this  is  a  measure  of  the 
mixed-layer  depth,  since  it  refers  to  the  thin  temperature  jump  zone  at  the 
base  of  the  mixed  layer.  The  depths  were  estimated  from  smoothed  contours 
of  the  thermistor-chain  data  such  that  any  variability  above  about  the 
inertial  frequency  does  not  appear.  Variability  at  and  above  the  inertial 
frequency  will  be  discussed  later. 

There  are  several  striking  features  about  Fig.  21.  Firstly,  as  mentioned 
in  the  previous  section,  despite  the  passage  of  four  storms  all  with  wind 
speeds  exceeding  10  m/s,  the  maximum  change  in  layer  depth  at  any  mooring 
was  still  only  about  half  that  of  the  south-north  (F-A)  change.  This  is  a 
measure  of  the  high  spatial  variability.  Secondly,  the  rate  of  deepening 
varies  considerably  with  mooring  position.  There  was  virtually  no 

deepening  at  position  D,  for  example,  while  during  the  same  period  at  E  the 
layer  deepened  by  about  20  m.  Such  a  result  is  a  certain  indication  of 
upwelling  at  or  upstream  from  position  D.  Thirdly,  some  considerable 
large-scale  variations  in  the  layer  depth  can  be  seen,  particularly  at 
mooring  position  F.  Here  between  14  and  17  Oct  the  depth  varied  from  63  m 
to  39  m  and  back  to  71  m,  with  a  similar  variation  occuring  later  between 
28  and  29  Oct.  Closer  examination  of  these  events  (to  appear  later) 

reveals  that  the  mixed  layer  itself  varied  by  somewhat  less  than  these 
extreme  values  and  that  the  variations  were  probably  due  to  advection 

effects  —  either  southward  advection  of  the  colder,  more  shallow  mixed 

layer  to  the  north,  or,  southeastward  advection  of  mixed  water  formed  in  a 
region  of  convergence  between  cold  northern  and  warmer  southern  waters. 

Figures  22  to  27*  present  a  detailed  day-by-day  description  of  the  tem¬ 
perature  structure  within  the  mixed-layer/thermocl i ne  zone  at  each  of  the 
six  mooring  positions.  The  contours  were  estimated  using  a  cubic  spline 
interpolation  (with  tension  factor  5)  between  thermistors  for  each  time 
step  and  then  using  a  linear  interpolation  in  time  between  successive  pro¬ 
files.  Upper  thermistor  and  lower  thermistor  temperature  values  are  writ¬ 
ten  at  their  corresponding  depth/time  positions  at  4  h  intervals  and,  since 
the  interval  between  isotherms  is  fixed  at  0.5°C,  the  absolute  value  of  an 
isotherm  is  readily  deduced. 

The  Nyquist  frequency  for  these  data  is  3  cycle/h,  but  since  the  declared 
time  response  for  the  thermistor  is  1  min  some  aliasing  has  clearly 
occurred.  However,  it  will  be  demonstrated  in  Sect.  5.1,  on  temperature 
spectra,  that  the  extent  of  this  aliasing  is  quite  limited  and  can  not 
significantly  affect  the  interpretation  of  these  contours. 

Various  periods  of  quite  strong  inertial  internal  waves  (of  18  h  period  at 
this  latitude)  can  be  seen  in  most  of  the  records.  These  vertical  oscilla- 


*The  to  tn  tht&  ahaptzA  OAZ  on  pp.  51 -SI. 
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FIG.  21  DEPTH  OF  MAXIMUM  THERMOCLINE  ESTIMATED  FROM  THERMISTOR  CHAIN  DATA. 
6-31  Oct,  1980. 
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tions  of  the  thermocline  are  quite  common  in  the  Mediterranean  and  usually 
follow  sudden  increases  in  wind  speed  (see  <14>  <15>  <16>,  for  example). 
Some  examples  are  evident  in  Fig.  24  for  mooring  C  on  11,  15  and  19  Oct. 
In  some  cases  inertial  oscillations  appear  to  follow  sharp  decreases  in 
wind  speed  —  for  example  on  22,  25  and  28  Oct  at  mooring  B,  as  shown  in 
Fig.  23.  In  other  cases,  however,  inertial  oscillations  seem  to  appear 
unforced  by  sudden  changes  in  wind  speed.  Two  such  examples  are  at  mooring 
E  (Fig.  26)  on  13  Oct  and  at  mooring  F  on  23  Oct  (Fig.  27).  It  is  possible 
that  the  local  winds  at  E  &  F  on  these  particular  days  did  indeed  change 

rapidly.  If  this  did  happen,  however,  one  would  expect  some  evidence  of 
the  variation  to  show  up  at  position  A,  only  150  km  to  the  north. 

We  will  see  in  Sect.  5.1  that  the  inertial  oscillations  are  well 

demonstrated  in  the  temperature  variance  spectra.  We  will  also  see  evi¬ 
dence  there  of  high  variance  at  the  unexpected  frequency  of  about  8  h. 
Since  this  does  not  correspond  to  either  the  inertial  or  tidal  frequencies, 
it  is  probably  due  to  some  non-linear  interaction  between  the  two.  A  good 
example  of  this  variance  can  be  seen  at  mooring  D  (Fig.  25)  on  14  Oct  which 
shows  a  fairly  energetic  oscillation  with  a  period  between  about  7.5  and 
8  h. 

Re-examining  the  large  vertical  excursion  of  the  thermocline  at  position  F 
shown  in  Fig.  21,  we  notice  that  during  14-17  Oct  (Fig.  27)  that  the  mixed 

layer  shallowed,  due  not  so  much  to  an  upward  movement  of  the  whole  ther¬ 

mocline  but  more  to  a  rarefaction  of  the  isotherms.  This  event  seems  most 
likely  to  be  the  result  of  southeastward  advection  of  waters  mixed  horizon¬ 
tally  somewhere  upstream  from  F  by  convergence  of  the  cold  northern  and 
warm  southern  waters.  A  likely  location  for  this  is  near  mooring 

position  E.  Assuming  a  cyclonic  eddy  to  the  north  and  anticyclonic  eddy  to 
the  south,  as  described  in  Sect.  2.2,  then  a  convergence  zone  would  be 
expected  near  position  E.  Indeed,  we  see  evidence  of  such  a  convergence 
during  14-16  Oct  (Fig.  25)  where  warmer  water  appears  to  be  advecting  in 
over  the  colder  water.  The  vertical  temperature  sections  taken  on  14  and 
15  Oct  (Fig.  8)  also  indicate  a  convergence  near  position  E. 


3.2  Current-meter  data 

Some  of  the  current-meter  data  were  presented  in  Sect.  2.4  in  the  form  of 
progressive  vector  diagrams  (Figs.  19,  20).  These  data  are  now  re¬ 

presented,  together  with  other  current-meter  data,  in  the  form  of  stick 
diagrams  in  Fig.  28  to  37.  The  vector  components  were  first  low-pass 
filtered  to  about  6  h  and  then  plotted  as  one  stick  every  hour.  This 
method  was  chosen  as  it  gives  a  better  impression  of  magnitude,  direction, 
and  sense  of  rotation  than  other  representations. 

Inertial  oscillations  can  be  seen  for  all  the  current  meters,  with  the 
strongest  occurring  nearest  the  surface.  Compare  Figs.  28  to  33,  which 
show  currents  measured  generally  below  the  thermocline  at  each  of  the 
mooring  positions,  with  Figs.  34  to  37,  which  show  currents  measured 
generally  within  the  mixed  layer  at  position  C.  The  current  measured  by 
the  NBA  current  meter  at  54  m  at  position  C  (Fig.  33)  appears  to  be  higher 
than  that  measured  at  32  m  by  the  VACM  current  meter  (Fig.  34).  This  is 
probably  due  to  the  pumping  of  the  NBA  current  meter  discussed  in 
Sect.  2.5. 
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A  Niel -Brown  3-axis  acoustic  current  meter  was  moored  at  4  m  depth  at 
position  C  in  an  attempt  not  only  to  measure  the  mean  current  at  this  depth 
(Fig.  44)  but  also  to  test  for  the  presence  of  Langmuir  circulations. 


3.3  Test  for  Langmuir  Circulations 

Langmuir  circulations  refer  to  the  vortex  structure,  sometimes  occurring  in 
the  mixed  layer,  of  counter-rotating  cell  pairs  with  axes  that  are  horizon¬ 
tal  and  aligned  approximately  in  the  wind  direction.  They  often,  but  not 
necessarily  always,  give  rise  to  the  parallel  bands  of  visible  substances 
collecting  in  the  convergence  lines  that  are  frequently  referred  to  as 
windrows.  The  most  commonly  accepted  theory  for  the  generation  of  Langmuir 
circulations  is  an  interaction  between  surface  waves  and  the  wind  drift 
current,  (see  <17>  <18>  <19>  <20>,  for  example).  Their  importance  is  in 
the  high  efficiency  with  which  they  transport  heat,  momentum,  and  turbulent 
kinetic  energy  downwards  through  the  mixed  layer  <18>,  but  little  is  known 
as  to  how  often  they  occur.  If  they  are  as  common  as  the  above-mentioned 
theories  seem  to  imply  they  would  consitute  a  serious  omission  in  the  pre¬ 
sent  state  of  mixed-layer  modelling.  On  the  other  hand,  should  they  occur 
only  in  certain  locations  and  on  the  rare  occasions  when  windrows  have 

actually  been  observed,  they  may  be  treated  as  anomalies  or  even  ignored  in 
their  over-all  effect  on  surface-layer  dynamics. 

Had  Langmuir  circulations  been  present  in  the  MILEX-80  zone,  their  vertical 
velocities  (up  to  3  or  4  cm/s  for  wel  1 -developed  cells)  ought  to  have  shown 
up  in  the  4m  vertical  velocity  measurement  of  the  Neil-Brown  current 
meter.  The  three  components  were  high-pass  filtered  to  3  h  and  rotary  and 
kinetic  energy  spectra  examined  for  peaks  in  the  2  min  to  3  h  range.  An 
advecting  Langmuir  circulation  structure  would,  at  some  frequency  related 
to  the  advection  speed  and  cell  dimensions,  cause  a  signal  in  the  rotary 

spectra  for  the  velocity  component  pairs  (u,v)  or  (v,w).  No  clear  signals 
were  observed  in  the  spectra  and  so  it  is  concluded  that  Langmuir  cir¬ 
culations  did  not  occur,  at  least  at  mooring  position  C,  during  the  11  day 

recording  period  of  the  Niel -Brown  current  meter. 


3.4  Wind  data 

Although  wind  speed  plots  have  been  included  in  many  of  the  figures  as  an 
indication  of  weather  conditions,  the  vector  wind  is  best  represented  by 
the  stick  diagrams  shown  in  Fig.  38.  This  shows  that  the  wind  was  blowing 
quite  consistently  from  the  east  throughout  the  experiment.  During 
7-13  Oct  the  wind  was  measured  on-board  ship  at  20  m  height  and  for  the 
remainder  of  the  experiment,  13-30  Oct,  the  wind  was  measured  at  a  height 
of  3.5  m  by  instruments  mounted  on  the  surface  buoy  at  A.  The  higher  value 
of  the  ship-measured  wind  is  due  partly  to  the  greater  height  of  the  ship's 
anemometer.  The  main  reason,  however,  is  that  mooring  A  was  more  sheltered 
than  the  zone  downwind  of  the  Bonifacio  Strait  where  the  ship  was  working 
during  7-13  Oct.  Indeed,  these  data  can  be  taken  as  further  evidence  that 
the  wind  is  much  stronger  off  the  Bonifacio  Strait  than  it  is  to  the  north 
and  it  is  this  variation  in  wind  stress  that  will  be  invoked  later  to 
explain  the  existence  of  the  cold  upwelled  water  in  the  north  of  the  zone. 
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3.5  Dynamic  stability 

Turbulent  entrainment  of  the  colder  water  in  the  thermocline  into  the  mixed 
layer,  and  hence  mixed-layer  deepening,  takes  place  if  the  local  Richardson 
number 


(Eq.  3) 


drops  below  0.25.  Here  a  is  the  coefficient  of  thermal  expansion,  u  the 
horizontal  velocity  vector,  T  the  temperature,  g  the  acceleratioF  due 
to  gravity,  and  z  the  vertical  axis.  This  mechanism  is  important,  for 
example,  in  the  Pollard,  Rhines  and  Thompson  model  <11>,  where  it  is 

assumed  that  the  mixed  layer  moves  as  a  slab  over  the  thermocline,  which  it 

then  entrains  when  the  shear  at  the  base  of  the  mixed  layer  becomes  suf¬ 
ficiently  large. 

Eddy  viscosity  models  of  mixed-layer  deepening  parameterize  the  turbulent 
diffusion  by  use  of  an  empirical  function  K  that  depends  on  R-j  and,  if 

higher  order  turbulence  closures  are  used  (<21>,  for  example),  on  q  ,  the 

turbulent  kinetic  energy.  One  very  simple  empirical  form  for  K  due  to  Munk 
&  Anderson  <22>  will  be  used  in  Sect.  4.5.  The  purpose  of  this  section  is 
to  present  estimates  of  the  Richardson  number  at  different  depths 
throughout  the  MILEX-80  experiment  so  that  periods  of  mixing  can  be  iden¬ 
tified.  These  estimates  were  made  using  data  from  the  five  current  meters 
and  thermistor  chain  at  position  C  and  then  using  Eq.  3  in  the  form 


(Eq.  4) 


where  N  is  the  Brunt-VaTsal a  frequency  and  S  the  shear.  The  values  of 
and  were  calculated  by  the  expressions 


(Ui  -  U2)^  +  (Vi  -  V2)^ 


(Eq.  5) 


(Zl  -  Z2)2 


with  the  suffixes  i  and  2  referring  to  the  upper  and  lower  current  meters 
respectively,  and 


N2  =  9  “(^n  ~  ^n+i) 


(Eq.  6) 


n  -  Zn+i 
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with  n  and  n+1  referring  to  the  two  thermistors  nearest  the  mid  point 
{zi+Z2)/2  between  current  meters.  Since  thermistors  were  closer  together 
than  current  meters  (see  Sect.  1.2)  the  value  of  S,  as  a  vertical  mean 
value  over  a  relatively  large  interval,  compared  with  N,  as  a  vertical 
mean  value  over  a  relatively  small  interval,  tends  to  cause  Eq.  3  to  be  an 
overestimate  of  Ri .  In  the  estimates  of  Rj  to  be  presented  in  this  sec¬ 
tion,  therefore,  turbulent  mixing  might  well  be  occurring  for  values  of  Rj 
greater  than  0.25.  For  example,  entrainment  is  defined  to  occur  in  the 
Pollard,  Rhines  &  Thompson  model  <11>  when  the  bulk  Richarson  number  falls 
below  unity.  Since  the  estimates  to  be  presented  here  lie  somewhere  bet¬ 
ween  the  bulk  and  local  Richardson  number,  the  choice  of  critical 
Richardson  number  ought  to  lie  somewhere  between  1.0  and  0.25. 

Figures  39  to  42  illustrate  the  estimates  of  Richardson  number  for  each  of 
four  depths,  together  with  the  wind  speed,  mixed-layer  depth,  estimates  of 
current  shear,  and  estimates  of  the  Brunt-VaTsal a  frequency  throughout  the 
MILEX-80  period.  Initially  the  mixed  layer  was  at  about  the  14  m  level 
(Fig.  39)  so  that  the  4  m  and  14  m  current  meters  were  both  within  the 

mixed  layer.  This  implies,  as  can  be  seen  in  the  plot  of  N,  that  the  stra¬ 
tification  was  low  but  that  the  shear  near  the  layer  base  was  high  and,  as 

expected,  the  Richardson  number  was  well  below  the  critical  value. 

Occasionally,  however,  the  mixed  layer  shallowed,  the  thermistor  at  14  m 
fell  within  the  thermocline,  and  the  values  of  N  and  R-j  both  increased, 
causing  greater  stability  at  this  level.  After  14  Oct  the  mixed  layer 

deepened,  the  shear  maintained  its  relatively  high  value  of  0.01  to 
0.02  rad/s  and,  apart  from  a  short  period  on  15  Oct,  the  4  to  14  m  layer 
became  one  of  turbulent  production. 

Figure  40  shows  N,  S  and  Rj  estimates  for  the  14  to  23  m  level.  During 

9  Oct  the  mixed  layer  still  had  not  deepened  beyond  18  m  and  the  estimates 
of  N  and  S  were  both  high,  indicating  the  high  temperature  and  velocity 
gradients  across  the  mixed-layer/thermocl ine  interface  zone.  It  is 
interesting  to  note  that  whenever  the  mixed  layer  dropped  below  the  lower¬ 
most  current  meter  the  value  of  S  decreased  to  its  much  lower,  mixed- 
layer  value  of  between  0.005  &  0.01  rad/s.  An  example  of  this  is  when  the 
mixed  layer  deepened  to  below  the  23  m  current  meter  on  18  Oct. 
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FIG.  22  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  A. 

Printed  numbers  indicate  temperature  and  position  for  the  upper 
and  lower  thermistors  at  4  h  intervals . 
a)  13-15  Oct  1980. 
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FIG.  22  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  A. 

Printed  numbers  indicate  temperature  and  position  for  the  upper 
and  lower  thermistors  at  4  h  intervals, 
h)  16-18  Oct  1980. 
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22  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  A. 
Printed  numbers  indicate  temperature  and  position  for  the 
and  lower  thermistors  at  4  h  intervals, 
c)  19-20  Oct  1980. 
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FIG.  23  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  B. 

Printed  numbers  indicate  temperature  and  position  for  the  upper 
and  lower  thermistors  at  4  h  intervals . 
a)  12-14  Oct  1980. 
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FIG.  23  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  B. 

Printed  numbers  indicate  temperature  and  position  for  the  upper 
and  lower  thermistors  at  4  h  intervals, 
b)  15-17  Oct  1980. 
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FIG.  23  BATHYTHERMS  MEASURED  BY  THE  THERMITOR  CHAIN  AT  POSITION  B. 

Printed  numbers  indicate  temperature  and  position  for  the  upper 
and  lower  thermistors  at  4  h  intervals . 
c)  18-20  Oct  1980. 
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FIG.  23  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  B. 

Printed  numbers  indicate  temperature  and  position  for  the  upper 
and  lower  thermistors  at  4  h  intervals, 
d)  21--23  Oct  1980. 
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FIG.  23  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  B. 

Printed  numbers  indicate  temperature  and  position  for  the  upper 
and  lower  thermistors  at  4  h  intervals . 
e)  24-26  Oct  1980. 
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23  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  B. 

Printed  numbers  indicate  temperature  and  position  for  the  upper 
and  lower  thermistors  at  4  h  intervals, 
f)  27 ’-29  Oct  1980. 
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FIG.  23  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  B. 

Printed  numbers  indicate  temperature  and  position  for  the  upper 
and  lower  thermistors  at  4  h  intervals . 
g)  30  Oct  1980. 
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FIG.  24  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  C. 

Printed  numbers  indicate  temperature  and  position  for  the  upper 
and  lower  thermistors  at  4  h  intervals . 
a)  7’-9  Oct  1980. 
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Fig.  24  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  C. 

Printed  numbers  indicate  temperature  and  position  for  the  upper 
and  lower  thermistors  at  4  h  intervals, 
b)  10-12  Oct  1980. 
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FIG.  24  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  C. 

Printed  numbers  indicate  temperature  and  position  for  the  upper 
and  lower  thermistors  at  4  h  intervals, 
c)  13-15  Oct  1980. 
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FIG.  24  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  C. 

Printed  numbers  indicate  temperature  and  position  for  the  upper 
and  lower  thermistors  at  4  h  intervals, 
d)  16-18  Oct  1980. 
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FIG.  24  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  C. 

Printed  numbers  indicate  temperature  and  position  for  the  upper 
and  lower  thermistors  at  4  h  intervals • 
e)  19-21  Oct  1980. 
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FIG.  24  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  C. 

Printed  numbers  indicate  temperature  and  position  for  the  upper 
and  lower  thermistors  at  4  h  intervals, 
f)  22-24  Oct  1980. 
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FIG.  24  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  C. 

Printed  numbers  indicate  temperature  and  position  for  the  upper 
and  lower  thermistors  at  4  h  intervals, 
g)  25-27  Oct  1980. 
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FIG.  24  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  C. 

Printed  numbers  indicate  temperature  and  position  for  the  upper 
and  lower  thermistors  at  4  h  intervals . 
h)  28  Oct  1980. 
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FIG.  25  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  D. 

Printed  numbers  indicate  temperature  and  position  for  upper  and 
lower  thermistors  at  4  h  intervals, 
a)  11-13  Oct  1980. 
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FIG.  25  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  D. 

Printed  numbers  indicate  temperature  and  position  for  upper  and 
lower  thermistors  at  4  h  intervals, 
b)  14-16  Oct  1980. 
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FIG.  25  MTHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  D. 

Printed  numbers  indicate  temperature  and  position  for  upper  and 
lower  thermistors  at  4  h  intervals, 
c)  17-19  Oct  1980. 
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FIG.  25  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  D. 

Printed  numbers  indicate  temperature  and  position  for  upper  and 
lower  thermistors  at  4  h  intervals . 
d)  20-22  Oct  1980. 
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FIG.  25  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  D. 

Printed  numbers  indicate  temperature  and  position  for  upper  and 
lower  thermistors  at  4  h  intervals, 
e)  23-25  Oct  1980. 
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FIG.  25  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  D. 

Printed  numbers  indicate  temperature  and  position  for  upper  and 
lower  thermistors  at  4  h  intervals, 
f)  26-28  Oct  1980. 
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FIG.  25  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  D. 

Printed  numbers  indicate  temperature  and  position  for  upper  and 
lower  thermistors  at  4  h  intervals, 
g)  29  Oct  1980. 
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FIG.  26  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  E. 

Printed  numbers  indicate  temperature  and  position  for  upper  and 
lower  thermistors  at  4  h  intervals . 
a)  10-12  Oct  1980. 
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FIG.  26  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  E. 

Printed  numbers  indicate  temperature  and  position  for  upper  and 
lower  thermistors  at  4  b  intervals . 
h)  13-15  Oct  1980. 
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FIG,  26  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  E. 

Printed  numbers  indicate  temperature  and  position  for  upper  and 
lower  thermistors  at  4  h  intervals, 
c)  16-18  Oct  1980. 
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26  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  E. 

Printed  numbers  indicate  temperature  and  position  for  upper  and 
lower  thermistors  at  4  h  intervals, 
d)  19-21  Oct  1980. 
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FJG.  26  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  E. 

Printed  numbers  indicate  temperature  and  position  for  upper  and 
lower  thermistors  at  4  h  intervals, 
e)  22-24  Oct  1980. 
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FIG.  26  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  E. 

Printed  numbers  indicate  temperature  and  position  for  upper  and 
lower  thermistors  at  4  h  intervals, 
f)  25-27  Oct  1980. 
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FIG.  26  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  E. 

Printed  numbers  indicate  temperature  and  position  for  upper  and 
lower  thermistors  at  4  h  intervals, 
g)  28  Oct  1980. 
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FIG.  27  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  F. 

Printed  numbers  indicate  temperature  and  position  for  upper  and 
lower  thermistors  at  4  h  intervals, 
a)  13-15  Oct  1980. 
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FIG.  27  BRTHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  F. 

Printed  numbers  indicate  temperature  and  position  for  upper  and 
lower  thermistors  at  4  h  intervals, 
b)  16-18  Oct  1980. 
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FIG,  27  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  F. 

Printed  numbers  indicate  temperature  and  position  for  upper  and 
lower  thermistors  at  4  h  intervals , 
c)  19-21  Oct  1980. 
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FIG.  27  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  F. 

Printed  numbers  indicate  temperature  and  position  for  upper  and 
lower  thermistors  at  4  h  intervals, 
d)  22-24  Oct  1980. 
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FIG.  27  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  F. 

Printed  numbers  indicate  temperature  and  position  for  upper  and 
lower  thermistors  at  4  h  intervals, 
e)  25-27  Oct  1980. 
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FIG.  27  BATHYTHERMS  MEASURED  BY  THE  THERMISTOR  CHAIN  AT  POSITION  F. 

Printed  numbers  Indicate  temperature  and  position  for  upper  and 
lower  thermistors  at  4  h  intervals, 
f)  28-29  Oct  1980. 
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SACLANTCEN  SR-75 


M1LEX80  NBA5S5  LPF6H 


FIG.  28  CURRENTS  MEASURED  BY  NBA565  AT  POSITION  B  (63  m) . 


MILEX80  NBA566  LPFGH 


FIG.  29  CURRENTS  MEASURED  BY  NBA565  AT  POSITION  D  (60  m). 
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SACLANTCEN  SR-75 


MILEX80  NBA5B7  LPF6H 


FIG.  30  CURRENTS  MEASURED  BY  NBA567  AT  POSITION  F  (88  m) . 


MILEX80  NBA5G8  LPFGH 


FIG,  31  CURRENTS  MEASURED  BI  NBA568  AT  POSITION  A  (58  m) . 
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SACLANTCEN  SR-75 


MILEX80  NBA570  LPFGH 


FIG.  32  CURRENTS  MEASURED  BY  NBA570  AT  POSITION  E  (54  m) . 


M1LEX80  NBA5G9  LPFGH 


FIG.  33  CURRENTS  MEASURED  BY  NBA569  AT  POSITION  C  (54  m). 
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SACLANTCEN  SR-75 


MILEX80  VACM357  LPFSH 


FIG.  34  CURRENTS  MEASURED  EiT  VACM357  AT  POSITION  C  (32  m) . 


MILEX80  VACM35B  LPFBH 


FIG.  35  CURRENTS  MEASURED  BY  VACM356  AT  POSITION  C  (23  m). 
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SACLANTCEN  SR-75 


MILEX80  VACM355  LPFBH 


FIG.  36  CURRENTS  MEASURED  BF  VACM355  AT  POSITION  C  (14  m) . 


MILEX80  NB3  LPFGH 


FIG.  37  CURRENTS  MEASURED  BY  NB3  AT  POSITION  C  (4  m) . 
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SACLANTCEN  SR-75 


MILEX80  MC707  WIND  ( LPF4H  ) 


FIG.  38  WINDS  MEASURED  El  METEO-CHAIN  MC707  AI  POSITION  A  (height  3.5  m) . 

Wind  from  7  to  13  Oct  1980  measured  by  ships  anemometer  at  20  m 
height . 
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FIG.  39  WIND  SPEED,  MIXED  LAYER  DEPTH,  CURRENT  SHEAR,  BUOYANCY  FREQUENCY, 
AND  RICHARDSON  NUMBER  AT  POSITION  C  (4  to  14  m). 

a)  8-13  Oct  1980. 

b)  14-20  Oct  1980. 
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FIG.  40  WIND  SPEED,  MIXED-LAYER  DEPTH,  CURRENT  SHEAR,  BUOYANCY 
AND  RICHARDSON  NUMBER  AT  POSITION  C  ( 14  to  23  m) » 
a)  7-12  Oct  1980, 
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FIG.  40  WIND  SPEED,  MIXED-LAYER  DEPTH,  CURRENT  SHEAR,  BUOYANCY 
AND  RICHARDSON  NUMBER  AT  POSITION  C  (14  to  23  m), 
b)  13-18  Oct  1980. 
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FJG.  40  WIND  SPEED,  MIXED-LAYER  DEPTH,  CURRENT  SHEAR,  BUOYANCY  FREQUENCY, 
AND  RICHARDSON  NUMBER  AT  POSITION  C  (14  to  23  m), 
c)  19-24  Oct  1980m 
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FIGm  40  WIND  SPEED,  MIXED-LAYER  DEPTH,  (CURRENT  SHEAR,  BUOYANCY  FREQUEN(CY , 
AND  RICHARDSON  NUMBER  AT  POSITION  C  (14  to  23  m), 
d)  25-30  Oct  1980. 
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FIG.  41  WIND  SPEED,  MIXED  LAYER  DEPTH,  CURRENT  SHEAR,  BUOYANCY  FREQUENCY, 
AND  RICHARDSON  NUMBER  AT  POSITION  C  (23  to  32  m) . 
a)  7-12  Oct  1980. 
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FIG.  41  WIND  SPEED,  MIXED-LAYER  DEPTH,  CURRENT  SHEAR,  BUOYANCY  FREQUENCY, 
AND  RICHARDSON  NUMBER  AT  POSITION  C  (23  to  32  m) . 
b)  13-18  Oct  1980. 
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FIG.  41  WIND  SPEED,  MIXED-LAYER  DEPTH,  CURRENT  SHEAR,  BUOYANCY  FREQUENCY , 
AND  RICHARDSON  NUMBER  AT  POSITION  C  (23  to  32  m) . 
c)  19-24  Oct  1980. 
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FIG.  41  WIND  SPEED,  MIXED-LAYER  DEPTH,  CURRENT  SHEAR,  BUOYANCY  FREQUENCY, 
AND  RICHARDSON  NUMBER  AT  POSITION  C  (23  to  32  m) . 
d)  25-30  Oct  1980. 
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FIG.  42  WIND  SPEED,  MIXED  LAYER  DEPTH,  CURRENT  SHEAR,  BUOYANCY  FREQUENCY, 
AND  RICHARDSON  NUMBER  AT  POSITION  C  (32  to  54  m) . 
a)  7-12  Oct  1980. 


FIG.  42  WIND  SPEED,  MIXED-LAYER  DEPTH,  CURRENT  SHEAR,  BUOYANCY  FREQUENCY , 
AND  RICHARDSON  NUMBER  AT  POSITION  C  (32  to  54  m) , 
b)  13-18  Oct  1980. 
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FIG.  42  WIND  SPEED,  MIXED-LAYER  DEPTH,  CURRENT  SHEAR,  BUOYANCY  FREQUENCY, 
AND  RICHARDSON  NUMBER  AT  POSITION  C  (32  to  54  m) . 
c)  19-24  Oct  1980. 
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FIG.  42  WIND  SPEED,  MIXED-LAYER  DEPTH,  CURRENT  SHEAR,  BUOYANCY  FREQUENCY, 
AND  RICHARDSON  NUMBER  AT  POSITION  C  (32  to  54  m) . 
d)  25-30  Oct  1980. 
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4  ATMOSPHERIC  FORCING 

4.1  Surface  heat  fluxes 


In  this  section  we  consider  changes  in  the  heat  content  of  the  surface 
layer.  We  will  show  how  the  estimated  surface  heat  fluxes  do  not  account 
for  the  observed  changes  in  heat  content  and  that  this  anomaly  must  be 
explained  by  the  occurrence  of  upwelling  and 'downwel ling.  An  Ekman  pumping 
mechanism  is  proposed. 

The  rate  of  heat  exchange  between  the  atmosphere  and  ocean  can  be  sub¬ 
divided  into  a  number  of  contributions  given  by 

0  =  Qs  +  Ql  +  Or  -  Qi  ,  (Eq.  7) 


where  the  subscripts  s,  1,  r  and  i  refer  to  sensible  heat  exchange, 
latent  heat  of  evaporation,  radiative  heat  loss,  and  irradiative  heating 
respectively.  The  first  contribution  to  Eq.  7  can  be  estimated  by  use  of 
the  bulk  aerodynamic  formula  (see  <23>  for  example) 


^s  "  ^T  Pa  ^p^’^s  "  ’’’a)  '^10  ’  (Eq.  8) 

where  U]^q  is  the  wind  speed  at  10  m,  Tg  is  the  air  temperature,  Tg 
the  sea  surface  temperature,  pg  the  air  density,  Cp  the  specific  heat  of 
air  at  constant  pressure,  and  Cj  an  empirically  determined  coefficient 
taken  here  to  be  0.002. 

Since  humidity  measurements  were  not  made  during  MILEX-80  the  evaporative 
heat  loss  Qj  of  Eq.  7  must  be  estimated  from  the  sensible  heat  loss  using 
the  Bowen  ratio  —  the  ratio  of  sensible  heat  fluxes  to  latent  heat  fluxes. 
The  Bowen  ratio  has  an  empirically  determined  form  dependent  mainly  on  the 
air/sea  temperature  difference  and,  to  a  lesser  extent,  on  the  wind  speed 
<24>).  For  the  range  of  interest  here  the  form  chosen  was 

0.16  -  0.089  (Tg-Ts),  [Tg-Tg  <-1.0] 

B  =  -  0.147(Tg-Ts)3  -  0.4(Tg-Ts)2  -  0.483(Ta-Ts ) ,  [-  KTg-Tg  <  0.5] 

0.16  -  0.89  (Tg-Ts).  [Ta-Ts  >  0-5] 

(Eq.  9) 


The  evaporative  heat  loss  is  then 

Ql  =  Qs/B.  (Eq.  10) 

The  non-radiati ve  heat  flux  Qs  +  Q-j  was  then  estimated  using  Eqs.  7,  8  X  9 
together  with  the  measurements  of  wind  speed,  air  temperature,  and  sea  sur- 
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face  temperature  recorded  by  the  meteo-chain  at  position  A.  It  was  assumed 
that  the  wind  speed  Uio  was  approximately  equal  to  the  measured  value  at 
3.5  m.  The  result  for  the  17  days  of  measurements  is  plotted  in  Fig.  43a. 
The  second  profile  of  this  figure  refers  to  the  non-radiati ve  heat  flux 
estimated  for  position  F.  This  estimate  is  based  on  measurements  of  wind 
speed  and  air  temperature  at  position  A  and  sea-surface  temperature 
measured  by  the  uppermost  thermistor  (37  m  depth)  at  position  F.  Although 
the  air  temperatures  at  A  and  F  should  be  approximately  the  same,  the  wind 
speed  at  F,  according  to  historical  data,  should  be  greater  than  at  A. 
Also,  the  37  m  temperature  is  a  deep  mixed-layer  value  that  will  generally 
be  less  than  the  surface  temperature.  These  underestimates  for  wind  speed 
and  sea-surface  temperature  therefore  imply  an  underestimate  for  the  rate 
of  cooling  at  F  and  the  profile  in  Fig.  43a  should  accordingly  be  regarded 
as  a  lower  bound. 

We  now  consider  the  radiative  heat  flux  contributions  Qp  and  Qj  in  Eq.  7. 
Unfortunately,  time  series  for  these  fluxes  cannot  be  accurately  estimated 
since  we  lack  a  detailed  knowledge  of  cloud  cover  during  the  experiment. 
However,  for  the  purposes  of  the  argument  to  be  developed  here,  it  will  be 
sufficient  to  show  that  spatial  variations  in  the  mean  (time  averaged) 
radiative  heat  fluxes  are  not  significant.  It  is  well  known  (see  <25>  for 
example)  that  climatological  averages  of  the  net  infrared  radiation  flux 
Qp  vary  strikingly  little  with  temperature.  We  do  not  expect,  therefore, 
that  the  south-north  gradient  in  sea-surface  temperature  in  the  MILEX-80 
zone  will  cause  a  significant  variation  in  the  radiative  heat  flux.  We 
shall  assume  a  spatial  uniformity  both  in  the  mean  cloud  cover  and  in  the 
surface  albedo  over  the  MILEX-80  zone. 

The  radiative  cooling  can  be  written,  after  Wyrtki  <26> ,  in  the  form 


Op  =  noT^  (0.39  -  0.05  (1  -  0.7  n)  +  4  naT|(Ts-Ta),  (Eq.  11) 


in  which  n  is  the  emissivity  of  sea  water,  usually  taken  to  equal 
0.985,  a  the  Stefan-Boltzmann  constant  equal  to  5.76  *  10“^  W/m^,  eg  the 
water  vapour  pressure  in  the  air  a  few  metres  above  the  sea  surface,  n  the 
fractional  cloud  cover,  and  Ts  and  Tg  the  temperatures  measured  in 
degrees  Kelvin.  The  term  in  T^  represents  a  small  correction  for 

stability  of  the  air  above  the  sea.  ^Since  the  humidity  was  not  measured  we 
must  again  use  the  Bowen  ratio  to  estimate  the  water  vapour  correction  term 
in  eg. 

By  approximations  involving  the  ratio  between  the  bulk  aerodynamic  for¬ 
mulae,  such  as  Eq.  8,  for  sensible  and  evaporative  heat  fluxes,  the  Bowen 
ratio  can  be  written  in  the  form  <27> : 


B  =  0.65 


(Ts-Ta) 

(es-ea) 


(Eq.  12) 


where  eg  is  the  water  vapour  pressure  in  the  air  at  the  sea  surface,  which 
is  assumed  always  to  be  saturated.  Using  an  integrated  form  of  the 
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NON-RAD!ATiVE  HEAT  FLUX 


RADIATIVE  HEAT  FLUX  COMPONENTS 


12  16  20  24  28  October  1980 

FIG.  43  COMPONENTS  IN  THE  SURFACE  HEAT  FLUX  FDR  POSITIONS  A  AND  F. 

(a)  non-radiative  (sensible  +  latent)  heat  flux,  (b)  coefficients 
in  the  expression  for  radiative  (incoming  and  out-going )  heat 
fluxes  (cloud  cover  corrections  not  included  -  see  text) 
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Clausius-Clapeyron  equation  (see  <28>  or  <29> )  eg  can  be  written  as  a  func¬ 
tion  of  temperature  as 


®s 


=  eo  exp 


(Eq.  13) 


The  values  Qq  and  Tq  of  Eq.  13  are  reference  values  found  experimentally  to 
be  6.11  mbar  and  273°K  respectively.  The  constant  L  is  the  latent  heat 
of  vapourization  for  water,  equal  to  2.55  *  10^  J/kg  at  20°C,  and  Ry  is  the 
gas  constant  for  water  vapour,  equal  to  461.5  J/kg/°K^ 

Equation  11  can  then  be  re-written  as 

Qr  =  qr  (1  -  0.7n2)  +  qc  ,  (Eq.  14) 

where 

qr  =  0.985  a  T^  (0-39  -0.5  e )  (Eq.  15) 

with  eg  calculated  by  use  of  Eqs.  13,  12  &  9  and 


qc  =  0.394  a  T3(Ts  -  Tg) 


(Eq.  16) 


Values  of  qr  and  qc  are  plotted  in  Fig.  43b  for  both  positions  A  and  F  and 

it  can  be  seen  that  in  comparison  with  the  non-radiati ve  heat  fluxes  they 

vary  insignificantly,  both  with  position  and  with  time. 

The  irradiance  from  the  sun  can  be  written  in  the  form  <30> : 

Qi  =  qi  (1  -  0.68  n)  ,  (Eq.  17) 

where  qi  is  a  daily  average  of  incoming  radiation  of  the  sea  surface  that 
depends  on  the  latitude  and  time  of  year.  Using  tables  obtained  by 

Berliand  <31>  and  discussed  in  Ivanoff  <32> ,  a  linear  form  for  qi  for  the 

month  of  October  and  at  41°N  can  be  extracted  and  expressed  by 


qi  =  237  -  1.88t 


(Eq.  18) 


with  t  the  time  in  days  of  October.  The  constant  varies  with  latitude 
by  about  5.4  W/m^  per  degree,  giving  rise  to  the  small  difference  observed 
in  qi  between  positions  A  and  F  plotted  in  Fig.  43b.  Time  integrals  of 
Eqs.  8,  10,  14  and  17  will  now  be  used  in  consideration  of  the  total  heat 
balance  for  the  surface  layer. 
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4*2  Surface-layer  heat  balance  and  upwellinq 

The  heat  balance  for  the  surface  layer  can  be  represented  by 

3  T  9  T 

P  p  Cp  i^.V^T  +  p  Cp  w  —  -  vV^T  =  0  ,  (Eq.  19) 


where  the  subscript  h  signifies  a  horizontal  vector,  z  is  positive 
upwards,  and  v  is  the  molecular  diffusion  coefficient.  Then,  replacing 

the  variables  w  and  T  by  w+w*  and  f+T'  respectively,  where  the  overbar 

represents  a  mean  over  several  hours  and  the  prime  a  turbulent  fluctuation, 
Eq.  19  can  be  averaged  and  integrated  vertically  to  give 

0  _  0  _  _ 

p  Cp  ^  /  T  dz  +  p  Cp  /  ^.Vh  T  dz 

z=-D  z=-D 


0  _  _ 

^  ^  ^  ^  ^  Cp(w'T')o  -  p  Cp(v^)_D 


0  _  0  _ 

-  P  Cp  /  Kh  v2  T  dz  -  p  Cp  V  /  v2  T  dz  =  0  ,  (Eq.  20) 

z=-D  z=-D 


in  the  derivation  of  which  the  continuity  condition 


^•ii  =  0  (Eq.  21) 

and  the  parameterization 


u'hT'  =  -khVhT  (Eq.  22) 

have  been  applied. 

The  first  term  of  Eq.  20  represents  the  rate  of  change  of  heat  content  per 
unit  area  for  the  surface  layer  of  depth  0,  the  second  term  the  horizontal 
advection,  the  third  term  the  vertical  advection  (upwelling  and 

downwelling),  the  fourth  term  the  turbulent  heat  flux  at  the  surface,  which 
is  set  equal  to  Q  of  Eq.  7,  the  fifth  term  the  turbulent  flux  of  heat 

through  the  lower  boundary  z  =  -D  ,  which  is  set  to  zero,  and  the  last  two 
terms  represent  the  horizontal  turbulent  diffusion  and  the  molecular  dif¬ 
fusion  of  heat,  both  of  which  can  be  neglected  in  comparison  with  other 

terms. 
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In  the  MILEX-80  data  we  shall  be  interested  in  the  mean  heat  balance  over  a 
13-day  period  and  so  we  can  introduce  a  long  time  average,  which  we  define, 
for  some  variable  (j),  as 


1  T 

<<!>>=  ±  /  4,dt  ,  (Eq.  23) 

t=0 

where  t  is  of  the  order  of  13  days.  The  balance  of  Eq.  20  then  reduces 
to 

A  H  ^  0  '  ' 

—  =  -  /  p  Cp  <Uh.VhT>  dz  -  <p  Cp  /  w  —  dz>  -  <Q>  ,  (Eq.  24) 

T  2=_d  ^  ^z=-D  3z 

where 

0  0 

AH  =  p  Cp  /  T(t)  dz  -  p  Cp  /  T(o)  dz 
z=-D  z=-D 


Assuming  that  the  cold  water  patch  is  caused  by  upwelling  and  not  by  hori¬ 
zontal  advection  from  some  distant  source,  the  measured  change  in  heat  con¬ 
tent  AH  well  within  the  cold  patch  will  be  independent  of  position. 
Without  loss  of  generality,  therefore,  we  can  choose  a  position  such  that 
the  mean  advection  is  zero  and  the  righthand  side  of  Eq.  24  reduces  to 
upwelling  and  surface-flux  terms.  If  now  we  consider  a  typical 

Mediterranean  mixed-layer  geometry  with  layer  of  depth  d,  a  region  of 
strong  thermocline  (often  referred  to  as  the  temperature  jump  zone)  of 
thickness  6  ,  and  then  a  region  of  progressively  weakening  temperature 
gradient  with  depth,  Eq.  24  can  be  written  approximately  as 


T 


“  <P  Cp  f 


■^-aj 

D  I  w  —  dz>  -  <Q> 
z=_d_(S  9z 


(Eq.  25) 


Defining  mean  representative  values  w*  and  r*  for  the  6-layer  such  that 

-d  _  -Y 

<  /  w  —  dz>  =  6  w*  r*  ,  (Eq.  26) 

z=-d-6 

we  can  write,  using  Eq.  25, 

^  -AH/t  -  <Q> 

P  Cp  6  r* 

The  mean  heat  flux  can  be  separated  into  two  components 

<Q>  =  <Qnr>  +  <Qr> 


(Eq.  27) 


(Eq.  28) 
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with  the  subscripts  NR  and  R  referring  to  non-radiati ve  and  radiative 
respectively.  Except  for  <Qr>,  all  the  parameters  of  Eq.  27  are  either 
known  or  can  be  estimated.  the  parameter  <Qr>  is  not  known  because  it  is 
strongly  dependent  on  cloud  cover,  which  was  not  measured  during  MILEX-80. 
However,  if  we  can  assume  that  <Qr>  does  not  vary  significantly  in  space 
then  we  can  make  an  estimate  of  the  difference  in  w*  between  two  locations. 
By  Eqs.  14  and  17  the  total  radiative  heat  fluxes  can  be  written 

Qr  =  qp  (1-0.7  n2)  +  qc  -  qi  (1-0.68  n)  .  (Eq.  29) 


The  coefficients  qp,  q^-  and  q-j  ^  plotted  in  Fig.  43b,  can  be  seen  to  vary  by 

insignificant  amounts  between  north  and  south  in  comparison  with  the  non- 
radiative  heat  flux  terms  (plotted  in  Fig.  43a).  Furthermore,  the  coef¬ 
ficients  qp  and  q-j  are  approximately  constant  in  time  for  the  13-day  period 
14-27  Oct.  The  difference  in  <Qr>  between  north  and  south  can  therefore  be 
reduced  to  the  expression 

<Qr>I\1  -  <Qr>s  =  0.7<qp>(<n2>s-<n2>fj)  -  0.68<qi>(<n>s-<n>N)  .  (Eq.  30) 


We  now  make  the  assumption  that  the  mean  and  variance  for  the  cloud  cover 
over  the  13-day  period  do  not  vary  significantly  in  space  and  the  righthand 
side  of  Eq.  30  can  be  equated  to  zero. 

The  difference  in  w*  between  north  and  south  can  therefore  be  written 


w 


* 


N 


(AHs-aHn)/t  +  <Q[\|r>$  -  <Qnr>N 
p  Cp  6  r* 


(Eq.  31) 


where  we  have  also  assumed  that  6  and  r*  have  not  varied  significantly 
between  north  and  south.  This  assumption  is  purely  for  convenience  since 
separate  values  for  r  and  6  for  north  and  south  are  available  from  the 
data.  All  the  quantities  on  the  righthand  side  of  Eq.  31  can  now  be  esti¬ 
mated. 


Mean  values  for  the  non-radi ati ve  heat  flux  are  estimated  by  integration  of 
the  profiles  in  Fig.  43a  over  the  13-day  period.  These  give 


<QnR>S  =  89  W/m2 
<QnR>N  =  53  W/m^ 


(Eq.  32) 


To  estimate  AH,  XBT  casts  were  divided  into  four  groups:  two  groups  in 
the  cold  water  patch  in  the  north  corresponding  to  the  beginning  and  end  of 
the  selected  13-day  period  and  two  similar  groups  in  the  warmer  water  to 
the  south.  The  temperature  profiles  were  integrated  from  a  depth  of  120  m 
to  the  surface,  giving  values  of 


0 

Hi  =  p  Cp  J  T-j  dz 
z=-120 


(Eq.  33) 
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which  were  then  averaged  within  each  group  to  give  an  estimate  of  surface- 
layer  heat  content.  Differences  over  the  period  for  the  north  and  south 
were  found  to  be 


AH^/t  =  -133  W/m2 
AHs/t  =  116  W/m2 


(Eq.  34) 


These  values  imply  a  net  gain  in  heat  in  the  south  that  could  be  due  to 
radiative  heating  or  downwelling  or  both. 

It  is  interesting  to  note  at  this  point  that  the  net  radiative  heat  flux  is 
negative,  meaning  that  insolation  exceeds  back  radiation  even  for  100% 
cloud  cover.  From  Eq.  29  and  values  deduced  from  Fig.  43b  the  mean 
radiative  heat  flux  for  the  relevant  period  can  be  written 

Or  =  -56  n2  +  136  n  -  100  (Eq.  35) 


For  n  =  1  (complete  cloud  cover)  Qr  =  -20  W/m^  and  for  n  =  0  (no  clouds) 
Qp  =  -100  W/m^.  Comparing  these  fluxes  with  the  non-radiative  fluxes  given 
in  Eq.  32  shows  that  there  is  a  range  of  mean  cloud  covers  that  would  pro¬ 
duce  a  mean  surface-heat  flux  of  zero. 

Using  the  values  of  Eqs.  32  and  34,  together  with  the  typical  values  of 
6  =  10  m,  r*  =  4°C/10  m,  q  =  1027  kg/m^ ,  and  Cp  =  1010  J/kg/°K,  the  mean 
vertical  velocity  difference  between  north  and  south  can  be  calculated  by 
Eq.  31  to  gi ve 

w*|\|  -  w*2  =5.9  m/day  (Eq.  36) 

Using  the  two  extreme  cases  Qr  =  -20  W/m^  and  Qp  =  -100  W/m^  deduced  from 
Eq.  35,  Eq.  27  shows  that  the  upwelling  velocity  w*f>j  in  the  north  lies  in 

the  range  2.1  to  3.8  m/day  and  the  downwelling  w*-  in  the  south  lies  in  the 
range  -3.9  to  -2.2  m/day.  ^ 

One  source  of  error  in  these  estimates  lies  in  the  use  of  the  Bowen  ratio 
technique  to  estimate  vertical  humidity  gradients  above  the  sea  surface. 
Other  error  sources  are  in  having  chosen  positions  A  and  F  for  the  estima¬ 
tion  of  the  heat  fluxes  (A  is  too  far  northwest  and  F  is  not  far  enough 
south),  in  the  assumptions  that  the  wind  speed  at  F  is  equal  to  the  wind 
speed  at  A,  and  in  the  use  of  a  deep  mixed-layer  temperature  to  approximate 
the  sea-surface  temperature  at  F.  By  far  the  greatest  source  of  error, 
however,  is  in  the  high  variance  of  the  estimates  for  heat -content  anomaly 
AH/t  made  from  XBT  profiles.  These  are  of  the  same  order  as  the  mean 
value  estimates  and  are  probably  due  both  to  insufficient  numbers  of  XBT’s 
used  in  the  estimates  and  to  poorly  defined  distinctions  between  north  and 
south  groupings. 
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4.3  Ekman  pumping 

In  this  section  we  propose  that  the  upwelling  and  downwelling  velocities 
estimated  in  the  previous  sections  are  due  to  Ekman  pumping  caused  by  the 
wind  jet  blowing  through  the  Strait  of  Bonifacio.  Historical  data  <10> 
show  that  for  most  of  the  year  the  wind  in  the  north  Tyrrhenian  Sea  blows 
from  the  west  and  that  the  intensity  of  this  wind  is  much  greater  in  the 
area  to  the  east  of  the  Bonifacio  Strait  than  it  is  either  to  the  north, 
which  is  partially  sheltered  by  Corsica,  or  to  the  south  which,  although 
similarly  sheltered  by  Sardinia,  appears  to  be  south  of  the  influence  of 
the  Mistral.  An  example  is  shown  in  Fig.  44.  The  figure,  taken  from  <10>, 
represents  the  summer  average  wind-stress  field  estimated  from  20  years  of 
shipping  wind  measurements. 

Figure  45  is  a  schematic  representation  of  the  Ekman  pumping  mechanism. 

When  the  wind  blows  over  the  sea  surface  the  top  layer,  in  which  the 

stresses  are  high  (in  our  case  confined  to  the  mixed  layer),  begins  to  move 
in  the  direction  of  the  wind.  Because  of  the  rotation  of  the  earth  the 
motion  deflects  to  the  right  and  after  a  time  this  results  in  a  net  mass 
flux  (or  Ekman  flux)  to  the  right  of  the  wind  direction.  When  the  wind 
varies  in  space,  as  it  does  across  the  jet  in  Fig.  45,  the  Ekman  flux  is 
greater  at  the  maximum  wind  stress  than  at  the  minimum  outside  the  jet. 

This  leads  to  a  divergence  zone  to  the  north  that  causes  upwelling  and  a 
convergence  zone  to  the  south  that  results  in  downwelling. 

Figure  45  also  shows  a  schematic  cross-section,  indicating  the  vorticity 

response  in  the  deeper  sub-surface  water.  This  response  will  be  invoked 
later  in  Sect.  4.4  to  explain,  in  part  at  least,  the  existence  of  the 
counter-rotating  meso-scale  eddies  observed  in  the  area.  Indeed,  it  is 
because  we  wish  to  examine  the  vorticity  imparted  to  the  water  by  the  wind 
that  we  derive  the  Ekman  pumping  expression  through  the  vorticity  balance 
rather  than  the  more  illuminating  description  through  the  momentum  balance. 

The  vertical  (normal)  component  of  the  vorticity  equation  on  the  rotating 
earth  can  be  written  in  non-dimensional  form  as  (see  App.  A  or  <33>  for 
further  explanation) 


Ro  ^  +  3v  +  Ro(ai..V)w  +  ^ 


Dt  -  8z 


_  1 


(Eq.  37) 
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SUMMER  WIND  STRESS  AVERAGE  FOR  POINTS  WITH  MORE  THAN  100  OBSERVATIONS 


SUMMER  - ^  1.0  DYNE/cm^ 


FIG.  44 


SUMMER  AVERAGED  WIND 


STRESS  OVER  THE  MEDITERRANEAN  < 10> . 


FIG.  45  SCHEMATIC  REPRESENTATION  OF  EKMAN  PUMPING  AND  SUBSURFACE  RESPONSE. 
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with  the  dimensionless  Rossby  number  Rg,  curvature  parameter  g,  and  hori¬ 
zontal  and  vertical  Ekman  numbers  Ej^  and  Ey,  respectively,  defined  by 


(Eq.  38) 


In  Eq.  37,  p  is  the  vertical  vorticity  relative  to  the  earth's  surface, 

-  oj  the  relative  vorticity  vector,  (x,y,z,)  a  right-handed  cartesian 
coordinate  system  with  z  the  upward  normal  axis,  (u,v,w)  tine 
corresponding  orthogonal  velocity  components,  p  the  density,  p  the 
pressure,  and  the  subscript  H  signifies  horizontal  components.  In  the 

relations  of  Eq.  38,  f  is  the  normal  component  of  the  earth's  vorticity 

(the  Coriolis  parameter),  gg  the  south-north  gradient  in  the  linear 
approximation  to  the  variation  of  f  with  latitude  (f  =  fg  +  ggy),  U  is 
a  horizontal  velocity  scale,  L  a  horizontal  length  scale,  D  a  vertical 
length  scale  (  «  L),  and  A^  and  A^  are  the  turbulent  eddy  viscosity 

coefficients  for  the  horizontal  and  vertical  respectively. 

The  fifth  term  of  Eq.  37  is  a  vorticity  production  term  due  to  the  barocli- 
nicity  of  the  fluid.  Although  this  term  is  formally  of  order  unity  it  has 
been  shown  by  Pedlosky  <33>  that  if  the  Rossby  number  is  small  (  «  1)  the 
value  of  this  term  will  also  be  small  and  can  be  neglected. 

In  our  case  the  horizontal  length  scale  is  that  associated  with  mean  south- 
north  variations  in  wind  stress,  which  is  about  100  km,  a  typical  velocity 

scale  for  the  Ekman  layer  is  about  10  cm/s,  and  the  local  value  of  f  is 

10“‘^/s.  These  give  a  value  of  Rg  =  0.01.  The  parameter  g  is  equal  to 
(cot  0g)L/rg,  with  6g  the  latitude  and  rg  the  earth's  radius,  and  has  the 
value  0.018.  Since  also  the  horizontal  Ekman  number  is  of  the  order  of 

D^/L^  less  than  the  vertical  (a  factor  of  1O“0  for  the  Ekman  layer)  we  are 

left  with  only  the  fourth  and  last  terms  of  Eq.  37. 

In  dimensional  form  and  integrating  over  the  depth  d  of  the  mixed  layer 
gives 


Rn  = 


g  = 


Eh  = 


fl 

3oL 

^0 

fL2 

^v 

fD2 


p  f  w 


z=-d 


,  9uh 

'H  ”  (->  ''v  z=o 


(Eq.  39) 


where  we  have  assumed  that  the  turbulent  stress  below  the  mixed  layer  is 
negligible  and  that  W2=g  is  zero.  The  term  in  brackets  is  the  horizontal 
turbulent-stress  vector  in  the  water  at  the  sea  surface  and  so  must  be 
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equal  to  the  tangential  wind  stress  t_.  The  upwelling  velocity  is  thus 
given  by 


w 


z=-d 


_  1 


pf 


curl  T 


Assuming  south-north  wind  variations  as  indicated  in  Fig.  45  this  reduces 
to 

1  3  T  Y 

w-d  =  -  ,  (Eq.  40) 

fp  8y 

where  y  is  the  northward  (or  cross  jet)  coordinate  and  x  the  eastward 
(or  downwind)  coordinate. 

The  mean  wind-stress  curl  over  the  northern  portion  can  be  written  as 


3y 

which  gives 


Pa  Cd 


(UlO+6Uio)^  -  (Uio-6Uio)^ 
L/2 


9y 


-  8  Pa  Cd  Uio  fiUio 


(Eq.  41) 


where  U^q  "le^P  wind  speed  over  the  region,  6Uio  the  wind  speed 

variation  amplitude  across  the  jet,  L  the  jet  wavelength,  pg  the  air 
density,  and  Cq  the  drag  coefficient.  Taking  Pa  =  1.22  kg/m^, 

Cq  =  1.5  *  10"3 ,  L  =  100  km.  Dig  =  8  m/s  and  6Uio  =  2  m/s  (i.e.  4  m/s 

difference  between  jet  maximum  and  minimum)  gives  "BT^/ay  =  -2.34*10"^N/m^ . 
The  mean  upwelling  velocity  is  then,  by  Eq.  40, 


w_c|  =  2.1  m/day 


(Eq.  42) 


This  value  is  to  be  compared  with  the  estimated  value  (range  2.1  to 
3.8  m/day)  of  the  previous  section.  In  view  of  the  historical  data  the 
requirement  of  a  wind  speed  4  m/s  higher  off  the  Bonifacio  Strait  than  off 
the  east  coasts  of  Corsica  and  Sardinia  appears  quite  possible. 


4.4  Sub-surface  response 

Descriptively,  the  sub-surface  layer  will  respond  to  Ekman  pumping  in  the 
following  way: 
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Upwelling  water  leaving  the  sub-surface  layer  will  cause  a  fluid  deficit 
(low  pressure)  and  hence  create  there  a  zone  of  horizontal  convergence. 
The  slowly  converging  water  will  deflect  to  the  right  of  its  motion  due  to 
the  earth's  rotation  and  thus  begin  to  circulate  in  a  cyclonic  motion 
around  the  upwelling  zone.  Similarly,  in  the  sub-surface  layer,  down- 
welling  water  will  cause  a  zone  of  horizontal  divergence  (under  high 
pressure)  with  the  outward  flowing  water  veering  constantly  to  the  right  of 
its  motion  and  circulating  in  an  anti  cyclonic  motion. 


To  what  extent  these  large  vortex  motions  can  grow  and  reach  the  propor¬ 
tions  observed  in  the  eddies  of  the  MILEX-80  zone  (i.e.  to  what  extent  the 
wind-stress  curl  is  a  generating  mechanism  for  these  eddies)  can  only  be 
determined  through  a  careful  examination  of  the  terms  in  the  vorticity 
balance,  Eq.  37.  For  simplicity  we  will  assume  a  homogeneous  fluid  and 
thus  avoid  the  problem  here  of  baroclinic  production  of  vorticity.  Also  we 
will  assume  a  flat  bottom  and  zero  vertical  velocity  at  the  top  and  bottomn 
boundaries.  If  at  the  lowest  order  a  geostrophic  balance  is  assumed  (i.e. 
the  Rossby  number  is  small),  then  as  shown  in  Pedlosky  <33>,  the  velocity 
components  u,  v,  and  hence  c  ,  will  be  independent  of  z  in  the  interior 
region.  Vertical  integration  of  Eq.  37  from  the  bottom  to  the  free  surface 
then  gives 


V 

Dt  Rq 


L_  +  Iv  3£ 
Ro  Rq_3z 


z=-l 


tqL 

pU2D 


VR  X  T_  , 


(Eq.  43) 


where  Tq  is  a  scale  surface  wind  stress  perturbation. 

The  third  term  of  Eq.  37  has  been  neglected,  since  it  can  be  written  as 


n  /su  3w  3v  3w 
*^0  [tz  ^  “  ■Jz  Tx 


+  4 


3w 

Jz 


and  since  u,  v  and  c  are  independent  of  z  to  the  lowest  order  in  the 
interior  and  w  is  zero  at  the  boundaries.  The  vertical  length  scale  D 
is  associated  now  with  the  total  depth  and  is  approximately  equal  to  1  km 
in  our  case.  We  shall  assume  that  the  time-averaged  velocity  near  the  bot¬ 
tom  is  in  reality  much  less  and  the  time-averaged  velocity  near  the  top  is 
in  reality  much  greater  than  the  average  velocity  over  the  whole  depth. 

Although  the  second  term  on  the  lefthand  side  of  Eq.  43  is  of  the  same 
order  as  the  first  it  will  be  small  if  spatially  averaged  over  the  eddy. 
The  principal  balance  in  this  case  would,  in  the  early  stages,  be  between 
the  tendency  for  the  vorticity  to  increase  (the  partial  derivative  contri¬ 
bution  to  the  first  term  on  the  lefthand  side  of  Eq.  43)  and  the  production 
by  the  curl  of  the  wind  stress.  Eventually,  of  course,  growth  would  be 
limited  by  friction  and  the  balance  would  be  between  the  wind-stress  term 
and  the  third  and  fourth  terms  on  the  lefthand  side  of  Eq.  43.  The  former 
balance  is  expressed  by  setting 

toL 

—5—  =  1  (Eq.  44) 

p  U^D 
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For  tq  =  0.025  N/m^  (a  perturbation  corresponding  to  2  m/s  wind  speed), 
p  =  1027  kg/m^ ,  L  =  100  km,  and  D  =  1  km  gives  a  velocity  scale  of 

U  =  0.05  m/s  .  (Eq.  45) 

Although  many  approximations  have  been  used  to  arrive  at  Eq.  45  and  any 
effects  of  bathymetry  and  mean  circulation  have  been  ignored,  this  value  is 
remarkably  high  considering  it  is  depth  averaged.  It  would  appear,  there¬ 
fore,  that  the  mean  wind-stress  curl  in  the  MILEX-80  zone  is  capable  of 
generating  the  two  counter-rotating  eddies  known  to  be  present  there. 
Furthermore,  generation  of  the  eddies  by  wind-stress  curl  would  account  for 
the  seemingly  high  barotropic  component  implied  by  the  current  measurements 
discussed  earlier  in  Sects.  2.4  and  2.5. 


4.5  Mixed-layer  modelling 

To  properly  model  the  mixed  layer  in  the  MILEX-80  zone  account  should  be 
taken  of  the  upwelling  velocities  due  to  the  Ekman  pumping  described  in 
Sect.  4.3.  A  model  representing  upwelling  effects  on  the  depth  of  the 
mixed  layer  was  developed  by  De  Szoeke  <34>.  The  model  assumed  very  large 
scale  variations  in  the  wind  stress,  as  found  in  the  Pacific  between  the 
Westerly  and  Trade  winds.  Here  we  can  apply  the  same  analysis  to  the  much 
smaller  scale  of  the  wind  jet  blowing  through  the  Bonifacio  Strait. 


We  consider  the  deepening  rate  of  the  mixed  layer  due  to  Kraus  and  Turner 
<6>  written  as 


3  /2 

dt  “  r  h^  a  g 


(Eq.  46) 


where  mg  is  an  empirically  determined  constant  (=1.2),  r  is  the  ther- 

mocline  gradient  below  the  jump  zone,  a  is  the  coefficient  of  thermal 
expansion  for  sea  water,  and  g  is  the  acceleration  due  to  gravity.  The 
mixed  layer  will  stabilize  at  some  depth  h_c|  where  the  rate  of  deepening 
(Eq.  46)  balances  the  Ekman  up-welling  given  by  Eq.  40.  Combining  these 
two  yields  ,  , 

\  / 


h_d 


4moT 


JagpV^  9Tx/3y 


(Eq.  47) 


Using  the  appropriate  values  T  =  0.5°C/m,  ag  =  0.002  m/s^/°C, 

T  =  0.1  N/m^  (i.e.  8  m/s  wind  speed),  and  3Tx/3y  as  in  Eq.  41  gives 

h_d  =  15  m  .  (Eq.  48) 

A  time  scale  for  the  adjustment  is  found  using  w_(j=  2.7  m/day,  estimated  by 
Eq.  42  to  give 

T  =  h_d/w_d  =  7  days  .  (Eq.  49) 
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The  value  of  15  m  is  close  to  the  measured  value  of  the  mixed-layer  depth 
found  in  the  cold  water  region  at  the  beginning  of  the  MILEX-80  experiment 
(end  of  summer  season)  and  the  time  scale  of  7  days  is  close  to  that  of 
passing  depression  systems  in  this  area.  Although  this  model  will  not  be 
developed  any  further  here  it  appears  to  be  the  right  sort  of  model  to 
explain  the  main  variations  in  the  depth  and  temperature  of  the  mixed  layer 
observed  in  the  north  Tyrrhenian  Sea. 


Other  areas  of  high  wind-stress  curl,  and  where  similarly  high  spatial 
variations  in  the  mixed  layer  might  therefore  be  expected,  can  be  seen  in  a 
plot  of  mean  annual  wind-stress  curl  due  to  May  <10>  shown  in  Fig.  46. 
These  areas  (and  associated  wind  systems)  are:  the  Gulf  of  Lions 

(Mistral);  northeastern  Agean  (Etesians);  western  Alboran  Sea  (Vendeval , 
Liberator  or  Levanter);  southern  Adriatic  or  northern  Aegean  (Sirocco, 
Gregale  and  Maestro);  northern  Adriatic  (Bora)  and  possibly  other  places 
not  showing  up  on  such  a  smoothed,  annually  averaged  and  non  synoptic  pic¬ 
ture.  Other  coastal  regions  not  in  the  Mediterranean  but  where  winds  are, 
as  in  the  Mediterranean,  channelled  onto  the  ocean  by  gaps  in  mountain 
ranges  may  also  be  areas  of  high  variability  in  the  mixed  layer. 

Despite  the  high  variability  in  the  MILEX-80  zone  shown  earlier  in 
Sect.  3.1,  both  in  the  plots  of  depth  of  maximum  thermocline  in  Fig.  21  and 
in  the  thermistor-chain  isotherms  plotted  in  Figs.  22  to  27  (variability 
caused  in  the  low-frequency  sub-inertial  range  by  advection  and  upwelling 
effects),  it  is  instructive  to  look  at  the  performance  of  a  one-dimensional 
mixed-layer  model. 


How  well,  for  example,  can  such  a  model  predict  the  general  trend  of  layer 
deepening  in  the  MILEX-80  zone  without  taking  into  account  the  effects  of 
upwelling  and  advection?. 


The  model 
empi 
dent 
model 


chosen  is  a  grid  point  or  differential  model  using  a  simple 
turbulence  closure  due  to  Munk  and  Anderson  <22>  that  is  depen- 


on 

the 

local  Richardson  number. 

The  diagnostic  equations  for 

the 

1  are: 

ly. 

-  fv 

a^u 

at 

az2 

+  fu 

a^v 

(Eq. 

50) 

at 

az2 

aT 

at 

=  Kt 

a^T 

az2 

9 

where  u  &  v  are  horizontal  velocity  components, 
dinate,  T  is  temperature,  and  K,,,  and  Kj 
ficients  for  momentum  and  heat  respectively. 


z  is  the  vertical  coor- 
are  eddy  viscosity  coef- 


The  turbulence  closure  is: 

K„  =  KJ(1  +  lORi)-'/" 
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Kj  =  K|(l  +  3.3Ri)'^/^ 


(Eq.  51) 


.  9T  , 

-a  g  —  / 

1  M 

9z 

9Z 

where  the  superscript  o  refers  to  a  constant  value,  R-j  is  the  local 

gradient  Richardson  number,  and  a  is  the  coefficient  for  thermal  expan¬ 
sion.  The  boundary  conditions  are,  at  the  surface: 


K  9  u 
'^m  — 

3z 


(Eq.  52) 


Kjll 

9z 


Qo 


and  below  the  thermocline: 

Km  If  =  (0.  0) 

.r  (Eq.  53) 

Kj  =  0  , 

9z 


where  and  xy  are  the  horizontal  components  of  the  surface  wind  stress 

and  Qp  is  the  surface  heat  flux.  The  estimates  of  wind  stress  and  heat 

flux  (non-radiati ve  component  only)  were  both  calculated  from  bulk  aerody¬ 
namic  formulae  (see  Sect.  4.1)  with  the  drag  coefficients  and  Cj  set  to 
0.002.  The  constants  and  Kj  were  both  set  to  0.02  m^/s.  Equations  49 

to  52  were  coded  numerically  using  a  "leap-frog"  numerical  scheme.  The 
model  was  run  for  a  17-day  period  and  compared  with  the  measurements  at 

position  C,  as  shown  in  Fig.  47. 

Agreement  between  the  model  output  and  the  general  trend  of  the  data  is 

quite  good  but  the  values  of  K,,,,  Kj,  Cfy,  and  Cj  required  for  this  agreement 

are  somewhat  high.  The  Munk  &  Anderson  turbulence  closure  <22>  gave  better 
results  than  another  form  tried,  since  it  allowed  for  higher  values  of  K,,, 

and  Kj  within  the  thermocline.  Looking  at  Fig.  28  it  can  be  seen  that  the 

reasonably  good  agreement  found  was  largely  due  to  the  choice  of  position  C 
for  the  comparison;  if,  for  example,  position  D  had  been  chosen  the 
agreement  would  have  been  poor. 

This  sort  of  model  might  be  used  successfully  in  an  area  such  as  MILEX-80 
if  it  were  adapted  to  include  an  upwelling  term  and  applied  at  a  number  of 
grid  points  in  the  horizontal.  Such  a  generalization  to  three  dimensions 
would  allow  the  upwelling  velocity  w  to  be  estimated  from  divergence  in 
the  Ekman  layer.  This  approach  has  been  tried  successfully  on  a  larger 
scale  by  De  Szoeke  <34>  using  a  bulk  mixed-layer  model. 
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CONTOURED  ANNUAL  WIND  STRESS  CURL 


FIG.  46  MEAN  ANNUAL  WIND-STRESS  CURL  OVER  MEDITERRANEAN  <10>. 


TCI 48  FILTERED  DATA  TEMPERATURE  (30h) 


12  14  16  18  20  22  24°C 


-«0- 


ML  DEPTH  TC148  (LPF30h] 


10  12  14  16  18  20  22  24  260ctober  80 


FIG.  47  COMPARISON  OF  1-D  MIXED-LAYER  MODEL  WITH  DATA  FOR  VERTICAL 
TEMPERATURE  PROFILES  (top  and  middle)  AND  FOR  MIXED-LAYER  DEPTH 
(bottom) . 
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5  SPECTRAL  ANALYSIS 


5.1  Thermistor-chain  spectra 

Before  presenting  the  temperature  variance  spectra  for  the  thermistor  chain 
data  it  is  appropriate  to  point  out  that  a  certain  degree  of  aliasing  is  to 
be  expected. 

Aliasing  can  occur  when  the  sampling  frequency  is  lower  than  the  frequency 
associated  with  the  time  response  of  the  measuring  instrument  (see  <35>, 
p.  228,  for  example).  The  frequencies  that  alias  (i.e.  fold  back  their 
energy  to  the  resolvable  spectrum  below  the  Nyquist  frequency)  are  given  by 
the  range 


1 

2At 


<  a 


(Eq.  54) 


where  At  is  the  sampling  frequency  and  6t  is  the  time  response  of  the 
instrument.  Since  the  specified  time  response  for  the  thermistor  is  1  min 
and  the  sampling  interval  was  set  at  10  min,  the  aliasing  range  is 
3  <  a  <  120  cycle/h. 

Examining  the  buoyancy  frequency  for  the  thermocline  at  position  C,  pre¬ 
sented  earlier  in  Sect.  3.5,  we  see  that  internal  waves  could  occur  in  the 
thermocline  at  frequencies  as  high  as  34  cycle/h  —  although  these  would 
have  low  energy.  As  the  aliasing  frequencies  lie  within  the  same  range  as 
the  high-frequency  internal  waves,  the  thermistor  frequency  response  was 
more  carefully  examined.  Departure  from  the  ideal  exponential  response  of 
the  thermistor  might  be  expected  to  be  caused  by  the  different  wrapping 
materials  used  in  forming  the  chain. 

An  impulse  of  finite  width  was  applied  to  the  thermistor  by  alternately 
plunging  the  instrument  into  two  baths  of  different  temperature  and 
measuring  the  impulse  response.  The  time  series  of  the  responses  were 
digitized  and  fourier  transformed  and  the  frequency  response  function  H  was 
calculated  from 


H(a)  =  I  I 
IGxlo)! 


(Eq.  55) 


where  Gxy(a)  is  the  cross  spectrum  of  the  impulse  x  with  the  response  y 
and  Gx  (a)  is  the  autospectrum  of  the  impulse. 

Figure  48  shows  the  impulse  and  response  time  series  for  three  different 
impulse  widths,  together  with  corresponding  frequency  response  functions 
estimated  according  to  Eq.  55.  The  high-frequency  cutoff  for  the  wrapped 
thermistor  appears  somewhat  better  than  the  120  cycle/h  value  assumed  from 
the  time  response.  At  10  cycle/h  the  frequency  response  is  only  0.25  and 
the  effective  aliasing  range  is  thus  reduced  to  between  3  and  10  cycle/h. 
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FIG.  48  THERMISTOR  FREQUENCY  RESPONSE  FUNCTIONS  (bottom)  ESTIMATED  FROM 
THREE  DIFFERENT  IMPULSES  (top). 
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Vertical  motions  of  the  mooring  line  occur  at  frequencies  well  above  this 
3  to  10  cycle/h  range  because  the  mooring  line  is  taut  nylon  and  the  buoy 
is  essentially  surface-following.  Internal  waves  in  this  frequency  range 
(and  buoyancy  frequency  range)  measured  by  Brekhovskikh  et  al  <36>  and  pre¬ 
sented  in  Phillips  <37>,  p.203),  show  vertical  displacement  spectra  of  the 
isotherms  decreasing  with  frequency  to  the  power  of  -2.7.  Similar  spectra 
for  a  lower  frequency  range  measured  by  Cairns  <38>,  also  presented  in 
<37>,  decrease  with  frequency  to  the  power  -2.  From  the  evidence  of  these 
two  spectra  any  aliasing  in  the  spectrum  from  the  internal  waves  in  the 

3  to  10  cycle/h  range  will  have  only  a  minor  effect. 

Temperature  variance  spectra  for  each  of  the  mooring  positions  A  to  F  are 
presented  in  Figs.  49  to  54. 

To  estimate  spectra  the  data  were  detrended  by  least  squares  and  multiplied 
by  a  cosine  taper  window  before  fourier  transforming.  Confidence  in  the 
raw  spectral  estimates  was  increased  by  averaging  the  raw  estimates  over 
adjacent  frequency  bands.  The  number  of  bands  was  set  to  seven  for  the 
spectra  of  the  individual  thermistors  and  to  three  for  the  averaged  spectra 
of  the  thermistor  chain.  Towards  higher  frequencies  the  number  of  bands 

averaged  was  progressively  increased  in  order  to  avoid  congestion  on  the 
logarithmic  frequency  axis.  Accordingly,  the  confidence  in  higher  fre¬ 
quency  estimates  is  also  increased,  as  can  be  seen  on  the  90%  confidence 
interval  limits  displayed  on  the  top  left  of  each  figure.  For  the  purpose 
of  display  most  of  the  spectra  have  been  cut  at  the  low-frequency  end  by 
linear  interpolation  between  the  first  two  estimates  (i.e.  all  spectra 
begin  at  the  frequency  of  0.01  cycle/h). 

Although,  as  noted  earlier,  aliasing  ought  not  to  be  serious  for  isotherm 
displacement  spectra,  it  does  appear  to  be  occurring  for  the  temperature 
variance  spectra.  It  appears  as  a  rapid  increase  in  variance  towards  the 
Nyquist  frequency,  as  clearly  evidenced,  for  example,  in  the  first  four 
thermistors  of  mooring  F  (Fig.  54).  Other  good  examples  occur  in  ther¬ 
mistor  T4  of  mooring  E  (Fig.  53)  and  in  thermistor  TIO  of  mooring  C 
(Fig.  50),  although  here  it  is  difficult  to  distinguish  between  aliased 

internal  waves  from  the  3  to  10  cycle/h  range  and  resolvable  internal  waves 
occurring  around  the  1  and  2  cycle/h  bands.  Aliasing  for  some  frequency, 
a  ,  is  caused  by  the  frequencies  2na[\j  +  a,  where  om  is  the  Nyquist  fre¬ 
quency  and  n  an  integer.  Taking  a  =  1  cycle/h  and  =  3  cycle/h  the 
aliasing  frequencies  will  be  5,  7,  11,  13,  ...  cycle/h.  Similarly  for 
a  =  2  cycle/h  the  aliasing  frequencies  will  be  4,  8,  10,  14  ...  cycle/h. 
Since,  as  clearly  shown  in  Fig.  49,  the  buoyancy  frequency  range  at  42  m 
depth  (depth  of  thermistor  TIO  for  mooring  C)  lies  within  the 

4  to  11  cycle/h  range,  and  since  the  peak  of  internal  wave  energy  is  at  the 

buoyancy  frequency,  we  deduce  that  at  least  one  of  the  two  peaks  (a  =  1, 
2  cycle/h)  in  the  TIO  spectrum  is  from  aliased  internal  waves.  It  is 
easily  seen,  however,  that  for  a  spectrum  proportional  to  the  aliased 
contribution  of  a  large  internal -wave  peak  at  about  6  cycle/h  would  be  vir¬ 
tually  indetectable  below  0.1  cycle/h.  In  support  of  this  all  the  spectra 

show  a  proportionality  of  0'',  with  r  between  -5/3  and  -2. 

The  averaged  spectrum  of  the  chain  is  much  smoother  than  the  separate  ther¬ 
mistor  spectra  and  is  therefore  a  good  indicator  of  the  spectral  peaks. 
However,  confidence  in  these  peaks  is  not  necessarily  increased  through  the 
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THERMISTOR  CHAIN  TEMPERATURE  VARIANCE 
SPECTRA  FOR  EACH  THERMISTOR  AND  FOR 
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FIG.  49 

THERMISTOR  CHAIN  TEMPERATURE  VARIANCE 
SPECTRA  FOR  EACH  THERMISTOR  AND  FOR 
CHAIN  AVERAGE  AT  POSITION  A. 
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FIG.  52 

THERMISTOR  CHAIN  TEMPERATURE  VARIANCE 
SPECTRA  FOR  EACH  THERMISTOR  AND  FOR 
CHAIN  AVERAGE  AT  POSITION  D. 


FIG.  51 

THERMISTOR  CHAIN  TEMPERATURE  VARIANCE 
SPECTRA  FOR  EACH  THERMISTOR  AND  FOR 
CHAIN  AVERAGE  AT  POSITION  C. 
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FIG.  54 

THERMISTOR  CHAIN  TEMPERATURE  VARIANCE 
SPECTRA  FOR  EACH  THERMISTOR  AND  FOR 
CHAIN  AVERAGE  AT  POSITION  F. 


FIG.  53 

THERMISTOR  CHAIN  TEMPERATURE  VTiRIANCE 
SPECTRA  FDR  EACH  THERMISTOR  AND  FOR 
CHAIN  AVERAGE  AT  POSITION  E. 
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averaging,  since  the  thermistors  of  the  chain  were  not  measuring  statisti¬ 
cally  independent  events. 

All  the  spectra  show  a  peak  at  the  inertial  frequency.  As  discussed  in 
Sect.  3.1,  strong  inertial  oscillations  occur  after  a  rapid  increase  in 
wind  speed.  These  have  been  observed  and  discussed  by  Pollard  <15>  and 
Gonella  <39>,  who  both  found  the  oscillations  to  be  confined  effectively  to 
the  surface  mixed  layer. 

A  theory  for  the  vertical  oscillations  at  the  inertial  or  near  inertial 
frequency  has  been  developed  by  Pollard  <15>.  In  the  models  explored, 
generation  of  the  vertical  oscillations  was  seen  to  depend  on  horizontal 
variations  (not  necessarily  local)  in  the  applied  wind  stress. 

In  the  following  we  emphasize  this  dependence  and  relate  the  amplitude  of 
the  vertical  oscillations  to  the  curl  of  the  local  wind  stress.  The  treat¬ 
ment  is  not  rigorous,  as  in  <15>,  but  is  presented  purely  to  illustrate  a 
mechanism  that  appears  to  be  operating  locally  n  the  MILEX-80  zone. 

In  the  thermocline,  temperature  fluctuations  at  the  inertial  frequency  will 
occur  primarily  through  advective  action  of  the  inertial  oscillations 
(tlH  ,  on  spatial  gradients  in  the  mean  temperature  field.  This  can  be 
represented  more  specifically  by  the  second  two  terms  of  the  temperature 
balance 


II 

at 


■‘■w-^  +  u^^v^  T  +  w  — 
8z  3z 


(Eq.  56) 


where  the  tilde  denotes  oscillation  at  the  inertial  frequency,  the  over-bar 
denotes  an  average  over  several  inertial  periods,  and  the  subscript  H 
refers  to  horizontal  components. 

Equation  56  can  now  be  used  to  relate  temperature  variance  at  the  inertial 
frequency  to  inertial  pulsations  in  the  Ekman  mass  flux  divergence  and, 
consequently,  to  the  local  wind-stress  curl.  To  develop  this  connection  it 
is  convenient  to  write  the  horizontal  advection  terms  of  Eq.  56  as 

UH  •  VhT  «  -  •  VO  77  ,  (Eq.  57) 

where  c  is  the  depth  of  a  sloping  isotherm  at,  say,  the  maximum  ther¬ 
mocline.  The  tilde  and  over-bar  symbols  will  be  re-introduced  later  but 
variables  in  Eq.  57  should  be  considered  low-pass  filtered  to  several 
hours. 

The  Ekman  mass  flux  divergence  is  obtained  by  vertical  integration  of  the 
continuity  equation 


•  iiH  +  =  0  (Eq.  58) 
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between  z  =  c  and  z  =  0  to  give 

Vh  •  Uh  =  -  _UH,^  •  »  (Eq.  59) 

where 

0 

Uh  =  /  UH  dz  .  (Eq.  60) 

z=c 


Combining  Eqs.  56,  57,  &  59  gives,  for  the  local  rate  of  change  of  tem¬ 
perature. 


II 

9t 


=  -  (^H  •  Uh)  ^ 


(Eq.  61) 


The  inertial  oscillations  satisfy  the  vertically  integrated  momentum 

bal ance 


^  Uh  +  f  k  X  jj|^  = 

at 


(Eq.  62) 


where  k  is  a  unit  vertical  vector,  the  horizontal  wind-stress  vector 
and  f  is  the  local  inertial  frequency.  The  extra  term  |^a?/at  from 

integration  of  the  time  derivative  in  Eq.  62  has  been  neglected,  since  it 
is  small  compared  with  the  terms  retained.  Taking  the  wind  stress  to  be  a 
step  function  applied  suddenly  at  time  t  =  0  and  to  be  in  the  x-direction 
only  (i.e.  west  wind  in  the  MILEX-80  zone),  leads  to  the  well-known  solu¬ 
tion 


U  =  —  sin  ft 
pf 

(Eq.  63) 

V  =  —  (cos  ft-1) 
pf 


that  defines  both  the  clockwise  rotating  inertial  oscillations  and  the 
Ekman  mass  flux.  Taking  the  divergence  of  Eq.  63  and  allowing  Xx  to  be  a 
function  of  y  only  (i.e.  we  assume  axx/ax  «  axx/ay)  gives,  for  the 
inertial  frequency  pulsations  in  the  Ekman  mass  flux  divergence. 


VH-ili 


1  ^''^x 

pf  ay 


cos  ft 


(Eq.  64) 
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Then  substituting  Eq.  64  into  Eq.  61  and  integrating  in  time  gives 


T  = 


1  ilx  3T 
f2p  8y  9Z 


sin  ft 


(Eq.  65) 


Finally,  equating  the  amplitude  term  in  brackets  with  the  measured  ampli¬ 
tude  from  the  spectrum  we  obtain 


il 

f2p  3y  8z 


1  /2 

2Aa  Sy(f) 


=  f2p  2Aa  Sj(f)  ^  /  Iy 


(Eq.  66) 


where  Sj(f)  is  the  spectral  density  estimate  of  temperature  variance  at 
frequency  f  and  Aa  is  the  frequency  bandwidth  of  the  estimate  before 
smoothing  in  frequency. 

Table  2  shows  the  modulus  values  of  wind-stress  curl  calculated  from  Eq.  56 
and  the  corresponding  Ekman  pumping  velocity  for  each  of  the  six  mooring 
positions.  The  value  Sj(f)  was  taken  as  the  chain  maximum  and  the  value 

3^-/3z  was  set  equal  to  the  representative  value  of  l°C/m.  The  wind-stress 
curl  is  higher  to  the  north  of  the  Bonifacio  Strait  than  to  the  south  and 
the  range  of  values  adds  strong  support  to  the  model  assumed  in  Sect.  4.3. 
The  values  in  the  table  do  not  distinguish  between  upwelling  and  down- 
welling  but  in  view  of  the  earlier  calculations  of  Sect.  4.2  it  seems  pro¬ 
bable  that  the  value  at  F,  at  least,  is  for  downwelling. 

It  should  be  pointed  out  that  the  solution  of  Eq.  66  represents  an  average 
of  the  wind-stress  curl  over  about  an  18-day  period  and  that  values  for 
short  periods  could  be  considerably  higher.  Also,  the  simple  model  chosen 
for  the  inertial  oscillations  was  for  convenience.  Other  models,  such  as  a 
slowly  growing  wind  stress,  could  equally  have  been  used. 

Before  leaving  temperature  variance  spectra  it  is  interesting  to  observe 
that  another  small  peak  is  evident  at  the  approximate  frequency  of  1/8 
cycle/h.  See,  for  example,  the  chain-averaged  spectra  for  mooring  posi¬ 
tions  B,  C  and  F  in  Figs.  50,  51,  and  54  respectively  or,  more  prominently, 
the  separate  thermistor  spectra  occurring  at  all  mooring  positions  but  par¬ 
ticularly  for  B  and  C  in  Figs.  50  and  51.  A  likely  explanation  for  peaking 
at  this  frequency  is  a  higher  order  mode  forced  by  non-linear  interaction 
either  through  inertial -inertial  or  tidal -inertial  components.  More  speci¬ 
fically  one  could  look  at  the  higher  order  solution  to  Eq.  61,  written  as 


^2 

3t 


(vh  .  Uh) 


111 

3Z 


(Eq.  67) 
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TABLE  2 

ESTIMATES  FROM  TEMPERATURE  VARIANCE  SPECTRA  FOR  MAGNITUDE  OF 
WIND-STRESS  CURL  1  |  AND  VERTICAL  EKMAN  PUMPING  VELOCITY  |co 


Mooring 

Positions 

A 

B 

C 

D 

E 

F 

depth  (m) 

26.5 

45.6 

32 

33.7 

56 

64.7 

s^(f)  •’cycpk 

15.6 

17.5 

15.9 

8.3 

6.1 

5.1 

sil  c  yn 

az 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

,  -3 

4cr  Cf>k^lo 

5.36 

2.33 

2.0 

2.28 

2.28 

2.65 

3.82 

2.67 

2.36 

1.82 

1.56 

1.54 

2.66 

1.86 

1.65 

1.27 

1.09 

1.08 
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where  Ti  -jg  defined  by  Eq.  65.  However,  this  investigation  will  not  be 
pursued  further  here. 

5.2  Current-meter  rotary  spectra 

Rotary  spectra  consist  of  clockwise  and  anticlockwise  components  giving  the 
sense  of  rotation  of  the  motion  at  each  frequency.  This  is  represented  in 
the  complex  plane  as  <39> 


Ua(t)  +  iv<j(t)  =  u-  e-Tot  +  u+  e+^ot  ^  58) 

a  a 


where  u^  and  v^  are  the  x  and  y  components  of  current  at  frequency 

a  and  u“  and  u'^  are  the  clockwise  and  anticlockwise  rotating 
a  a 

components  also  at  frequency  a  . 

Equating  sine  and  cosine  terms  in  the  fourier  representation  of  the 
lefthand  side  of  Eq.  68  allows  the  rotary  components  to  be  written  in  terms 
of  the  conventional  fourier  coefficients  as 


+ 

Ua 


Ua 


where  the  subscripts  u  and  v 
rotary  spectra  components  are  then 


_*  _ 

Ua 

~2Ar 


-  bu)(j 

^u)a  ,  (Eq. 

refer  to  the  velocity  components, 
given  by 


69) 

The 


+*  + 
~ZKF~ 


where  Aa  is  the  bandwidth  of  the  estimate.  The  rotary  coefficient  is 
defined  by 


Ra  = 


Ao  -  A^ 


a: 


^  a; 


(Eq.  71) 
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The  spectra  were  computed  and  smoothed  by  averaging  estimates  in  adjacent 
frequency  bands,  as  described  in  Sect.  5.1.  Figures  55,  56,  and  57 
illustrate  the  rotary  spectra  and  coefficients  for  all  the  MILEX-80  current 
meters.  The  spectrum  labelled  VACMSNB3  in  Fig.  57  is  the  average  spectrum 
over  all  VACM  and  NB3  current  meters.  The  one  labelled  NBA-567,  also  in 
Fig.  57,  is  the  average  over  all  NBA  current  meters  not  including  NBA567  at 
mooring  F.  All  spectra  show  a  peak  of  energy  at  the  inertial  frequency  in 
the  clockwise  component  and,  curiously,  for  the  NBA  current  meters  a  lower 
energy  but  still  a  significant  peak  in  the  anticlockwise  component. 
Whether  or  not  this  latter  peak  is  an  artifact  of  the  mooring  system  (NBA 
current  meter  on  taut  surface  mooring  —  see  Sect.  2.5),  which  seems  unli¬ 
kely,  can  only  be  tested  by  comparison  with  other  current  meters  and/or 
with  subsurface  moorings.  Such  a  comparison  cannot  be  made  effectively  for 
the  MILEX-80  data  since  all  moorings  were  of  the  taut-nylon  surface-float 
type  and  the  VACM  current  meters  were  generally  within  the  mixed  layer, 
whereas  the  NBA  current  meters  were  generally  below  the  thermocline. 

Another  possibility  follows  a  similar  argument  to  that  in  Sect.  5.1  con¬ 
cerning  the  inertial  peak  in  the  temperature  variance  spectra.  One  could 
consider  a  two-layer  system  with  the  upper  mixed  layer  diverging  horizon¬ 
tally  due  to  the  wind-stress  curl  and  a  lower  homogeneous  layer,  of  a 
higher  density,  responding  to  the  resultant  vertical  velocity  induced  at 
the  interface.  Such  a  two-layer  geometry  fits  quite  well  in  the 
Mediterranean,  especially  during  the  autumn  period  when  surface  wind  mixing 
greatly  sharpens  the  thermocline.  The  appropriate  momentum  balance  for  the 
lower  layer  is  then  given  by 

H  +  f  X  u  =  -  i  Vp  .  (Eq.  72) 

9t  p 


The  free  surface  responds  to  the  doming  of  the  interface  such  that  in  the 
steady  state  the  mean  horizontal  pressure  gradient  in  the  lower  layer  is 
zero.  The  dynamic  response  is  a  mean  cyclonic  motion  in  the  upper  layer. 
The  lower  layer  will,  however,  experience  temporal  fluctuations  in  pressure 
gradient  due  to  the  pulsating  Ekman  divergence  in  the  surface  layer.  An 
expression  for  the  pressure  at  depth  z  =  zq  is 

p  =  -  g  P2  ZD  +  (P2  -  Pi)  gc  .  (Eq.  73) 


where  P2  is  the  density  of  the  lower  layer,  pi  the  density  of  the  upper 
layer,  and  ^  the  z-coordinate  of  the  interface.  Setting  P2  -  p^  = 
gaAT,  where  a  is  the  coefficient  for  thermal  expansion  and  AT  the 
temperature  difference  between  layers,  the  righthand  side  of  Eq.  73  becomes 

-  1  Vp  =  -  (gaAT)  Vc  -  (gaV  Ti)c  ,  (Eq.  74) 

P 


where  Ti  is  the  mixed-layer  temperature  and  the  lower-layer  temperature 
is  taken  to  be  constant  in  space.  Then  by  use  of  the  kinematic  boundary 

condition 
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FIG.  55  ROTARY  SPECTRA  AND  COEFFICIENTS  FDR  NBA  CURRENT  METERS  AT  MOORINGS 
(from  top  to  bottom)  B,  D,  F,  AND  A. 
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FIG.  56  ROTARY  SPECTRA  AND  COEFFICIENTS  FOR  NBA  CURRENT  METERS  AT 
MOORINGS  C  (top)  AND  E,  THE  NIEL'  BROWN  3- AXIS  CURRENT  METER  AT 
MOORING  C  AND  THE  UPPER  MOST  VACM  CURRENT  METER  ALSO  AT  C 
(bottom) . 
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FIG.  57  FOTARY  SPECTRA  AND  COEFFICIENTS  FOR  VACM  CURRENT  METERS  AT 
MOORING  C,  THE  AVERAGE  FOR  VACM  AND  NIEL  BROWN  CURRENT  METERS  AT 
C  AND  (bottom)  THE  AVERAGE  FOR  ALL  NBA  CURRENT  METERS  EXCLUDING 
MOORING  F. 
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w?  -  UH,?-VC  =  ^ 


(Eq.  75) 


and  Eq.  59  we  can  write 


t  t 

_  1  Vp  =  -  (gaAT)  V  (/  Vh.^  dt)  -  gaVT  /  Vh-Uh  dt 


(Eq.  76) 


Substituting  for  the  Ekman  mass  flux  divergence  given  by  Eq.  64  and,  for 
simplicity,  assuming  that  both  the  mixed-layer  temperature  and  the  wind 
stress  are  functions  of  y  only,  we  have 


gg 

-  f2p 


9tx 
9y  9y 


+ 


S2L 

f2p 


AT 


9y2 


sin  ft 

y 


(Eq.  77) 


On  the  righthand  side  of  Eq.  77  the  ratio  of  the  first  term  to  the  second 
term  can  be  expressed  as 


9Ti  I 

R  =  -  .  — ,  (Eq.  78) 

9y  ZttAT 


where  L  is  the  wavelength  for  variations  in  the  wind  stress.  Taking 
9Ti/9y  »  1°C  in  20  km  (see  Fig.  11a),  L  =  100  km  and  AT  =  4°C  gives  a 
value  R  =  0.2.  Therefore,  neglecting  the  first  term  on  the  righthand  side 
of  Eq.  77,  Eq.  72  can  be  written 


9t 

—  +  f2  V  =  -  -2“  at  cos  ft  .  (Eq.  79) 

9t2  fp  9y^ 

The  equation  for  v  is  one  of  forced  undamped  oscillations,  which,  with 
the  initial  condition  u=v=0  att=0,  has  the  solution 

V  =  Bt  sin  ft 


u  =  -  Bt  cos  ft  +  —  sin  ft  (Eq.  80) 

f 
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g  _  ggAT 

2f2p  3y2 


In  terms  of  the  rotary  components  defined  by  Eq.  69  the  solution  is  written 
as 


u 


+ 

f 


B 

2f 


i 


u-  =  Bt  +  —  i 
f  2f 


(Eq.  81) 


The  term  in  t  tends  to  increase  clockwise  rotating, kinetic  energy  with 
time.  For  small  t  the  inertial  oscillations  are  not  rotating,  but  by 
t  -  l/2f  (1.5  h),  u^  =  2u;^  and  the  clockwise  energy  is  four  times  greater 

than  the  anticlockwise. 

Although  the  above  model  is  capable  of  distributing  inertial  energy  between 
clockwise  and  anticlockwise  components,  it  would  require  unrealistically 
high  values  of  the  second  spatial  derivative  in  wind  stress  if  it  were  the 
only  contribution  to  the  anticlockwise  peak  observed  in  the  rotary  spectra. 
For  example,  substituting  Eq.  81  into  Eq.  70  yields 


3y2 


4f3p 

-  X 

gaAT 


(2Aa  A+)^/^ 


(Eq.  82) 


which,  taking  the  values  Aa  =  0.007  cycle/h  and  A+  =  800  (cm/s)2/cycle/h 
from  the  spectrum  of  NBA565  at  position  B  (Fig.  55)  gives  the  value 


9 

- -  =  1.28  *  10-8 

3y2 


(Eq.  83) 


Using  the  estimate  of  3Tx/3y  =  2.6*10-6N/m3  estimated  from  the  temperature 
variance  spectrum  for  position  B  (see  Table  2)  and  the  value  given  in 
Eq.  83,  gives  a  length  scale  for  Tx  =  1.27  km.  Such  a  value  is  too  small 
by  a  factor  of  about  100. 

In  conclusion,  a  pulsating  Ekman  mass  flux  divergence  in  the  surface  layer 
is  not  sufficient  to  explain  the  observed  peaks  in  the  anticlockwise  com¬ 
ponents  of  the  rotary  spectra.  However,  this  might  be  a  consequence  of  the 
chosen  distribution  for  t  and  T][  (i.e.  a  dependence  on  y  only).  A 

dependence  of  these  variables  on  x  (the  downwind  direction)  might  intro¬ 
duce  an  anticlockwise  velocity  component  that  is  also  growing  in  time. 
Because  of  the  very  low  dissipation  below  the  thermocline,  resonant  forcing 
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of  the  type  defined  by  Eq.  79  could  cause  inertial  oscillations  (both 
clockwise  and  anticlockwise)  with  very  high  energy  to  develop  in  the  lower 
1 ayer. 


5.3  Wind  and  current  coherence  spectra 

Wind  speed  has  been  plotted  variously  in  earlier  sections  of  the  report 
(see  for  example  Fig.  40)  and  wind  vectors  were  represented  by  the  stick 
diagrams  shown  in  Fig.  38.  Similarly,  plots  of  current  are  given  in 
Figs.  19  and  20  in  the  form  of  progressive  vector  diagrams  and  in  Figs.  28 
to  37  in  the  form  of  vector  sticks.  In  this  section  we  look  at  wind/ 
current  coherence  and,  for  the  lowest  resolvable  frequency  band,  current/ 
current  coherence. 

When  estimating  coherence  functions  it  is  necessary,  in  order  to  obtain 
sufficient  confidence  in  the  estimate,  to  average  over  a  large  number  of 
adjacent  frequency  bands.  Alternatively,  the  averaging  could  be  over  the 
estimates  taken  from  a  large  number  of  separate  pieces  of  the  same  time 
series.  However,  where  high  energy  dominates  at  low  frequencies  (i.e.  when 
spectra  are  red)  the  lowest  frequency  estimate  for  each  piece,  and  hence 
also  for  the  average,  can  be  very  distorted  unless  a  high-pass  filter  with 
a  very  narrow  transition  bandwidth  were  applied  to  the  original  time  series 
before  division  into  pieces.  Since  the  wind  record  is  very  short  (only  18 
days)  and  the  storm  period  is  between  3  and  5  days,  any  filter  passing  the 
storm  frequency  would  require  so  many  weights  that  pieces  near  the  ends  of 
the  time  series  would  be  very  poorly  filtered.  Furthermore,  the  lowest 
frequency-estimate  when  averaging  over  pieces  contains  no  contribution  from 
frequencies  lower  than  Fp  =  1/Tp,  where  Tp  is  the  piece  length.  The  fre¬ 
quency  averaged  estimate,  on  the  other  hand,  represents  contributions  from 
frequencies  as  low  as  Ff  =  1/NTp,  where  N  is  the  number  of  pieces  in  the 

piece-averaged  estimate.  Although  these  low-frequency  contributions  to  the 
lowest  frequency  estimate  are  also  distorted  (frequencies  lower  than  Ff 
could  not  have  been  filtered  out)  frequency  averaging  was  adopted  as  the 
better  procedure  for  all  the  coherence  estimates  presented  in  this  section. 

Wind-speed/current-speed  coherence  and  phase  spectra  are  presented  in 
Figs.  58  and  59.  In  all  cases  the  wind  speed  is  that  measured  by  the 
meteo-chain  at  mooring  A  and  covers  the  period  13-30  Oct  1980.  Both  wind 
and  current  records  were  low-pass  filtered  and  decimated  in  time  to  15-min 
intervals.  Coherence  estimates  were  then  made  for  the  overlapping  portions 
of  each  record  pair  by  taking  the  ratio  of  the  frequency-averaged  cross 
spectrum  to  the  product  of  the  frequency-averaged  spectra.  Averaging  was 
over  a  minimum  of  twenty  adjacent  frequency  bands  in  all  cases  except  for 
estimates  involving  the  current  at  F  (NBA567,  Fig.  58  and  the  current  at 
C  (VACM355,  Fig.  59),  where  averaging  was  over  fifteen  adjacent  frequen¬ 
cies,  and  the  current  at  C  (NB3,  Fig.  59)  with  averaging  over  only  eight 
adjacent  frequencies.  The  averaging  was  chosen  to  give  the  highest 
possible  confidence  in  each  estimate  but  still  to  retain  in  each  case  one 
estimate  covering  the  range  below  but  not  including  the  inertial  frequency. 

All  the  spectra  show  a  low  coherence  in  the  second  (no  inertial  frequency 
components  in  the  wind)  and  higher  frequency  estimates.  The  highest 
coherence  in  each  case,  appart  from  those  with  current  meters  at 
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COHERENCE 

(WIND  SPEED, CURRENT  SPEED)  PHASE 

MC  707  (A), 


FREQUENCY  ( CYC/HR  ) 


FIG.  58  COHERENCE  AND  PHASE  SPECTRA  COMPARING  WIND  SPEED  AT  POSITION  A 
WITH  CURRENT  SPEED  AT  POSITIONS  B,  D,  F,  A  AND  C. 
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COHERENCE 

(WIND  SPEED, CURRENT  SPEED)  PHASE 

MC  707  (A), 


FIG.  59  COHERENCE  AND  PHASE  SPECTRA  COMPARING  WIND  SPEED  AT  POSITION  A 
WITH  SPEED  AT  POSITIONS  E  AND  C  (4  CURRENT  METERS  AT  C) . 


118 


SACLANTCEN  SR-75 


position  C,  is  for  the  lowest  frequency  estimate  that  represents  contribu¬ 
tions  from  the  approximate  range  from  19  h  to  17  days.  The  best  coherence 
is  found  between  the  wind  and  current  measured  at  position  A.  This  shows 
values  higher  than  0.5  for  periods  as  short  as  3  h.  Curiously,  the  worst 
coherence  is  found  at  position  C  (wind  at  A,  current  at  C),  even  for  the 
NBA  current  meter,  which  was  at  a  similar  depth  to  all  the  other  NBA 

current  meters.  Wind/current  coherence  is  virtually  insignificant  for  the 
NB3  and  the  three  VACM  current  meters  at  position  C  (Fig.  59)  and  perhaps 
it  is  of  significance  that  these  current  meters  were  within  the  mixed  layer 
for  most  of  the  experiment.  That  coherence  is  so  much  lower  at  C  than  at  A 
lends  support  to  the  assumption  that  the  wind  speed  off  the  Bonifacio 
Strait  is  different  from  that  to  the  north,  which  is  sheltered  by  Corsica. 
Small  peaks  in  the  wind/current  coherence  standing  out  against  the 
background  are  observed  at  1  cycle/h  both  for  the  NB3  current  meter  of 
position  C  moored  at  4  m  (Fig.  59)  and  to  a  lesser  extent  for  the  NBA 
current  meter  of  position  F  moored  at  88  m  depth  (Fig.  58).  Current  meter 
VACM356,  also  at  position  C,  shows  a  similar  peak  in  its  coherence  with 

wind  speed  at  0.79  cycle/h. 

The  generally  higher  coherence  at  low  frequency  can  be  attributed  to  the 
response  of  the  subsurface  current  to  the  storm  frequency  in  the  3  to  5  day 
range.  In  support  of  this,  an  unsmoothed  cross  spectrum,  not  shown  here, 
showed  evidence  of  peaking  in  the  same  3  to  5  day  range.  The  wind  clearly 
contained  very  low  variability  at  the  inertial  frequency,  evidenced  by  the 
very  low  coherence  seen  in  the  second  estimate  for  each  of  the  spectra. 

A  complication  in  interpreting  these  results  is  the  high  probability  of  the 
NBA  current  meters  being  pumped  by  surface  waves  (see  Sect.  2.5).  Since 

the  sea  state  also  responds  to  the  storm  frequency  it  is  not  possible  to 

separate  the  effects  of  surface  waves  and  subsurface  response  on  the  high 
coherence  observed.  It  is  unfortunate  that  all  of  the  deeper  current 
meters  were  NBA's. 

Estimates  were  also  made  of  the  coherence  in  current  speed  between  all 
pairs  of  current  meters  for  the  range  from  19  h  to  17  days.  The  values, 
including  those  with  wind  speed,  are  shown  in  Table  3.  The  values  with 
asterisks  represent  the  higher  frequency  range  of  13  to  25  h,  the  lower 
frequency  range  being  too  low  for  the  short  NB3  time  series.  For  this 
higher  range,  which  includes  the  inertial  frequency,  a  high  coherence  can 
be  seen  between  the  two  uppermost  current  meters  at  position  C  (4  m  & 
14  m).  Surprisingly,  however,  there  is  no  coherence  in  this  range  between 
the  4  m  and  23  m  or  the  4  m  and  32  m  current  meters.  For  this  same  range, 
but  not  shown  on  the  table,  there  is  a  coherence  of  0.51  between  the 
VACM355  (14  m)  and  VACM356  (23  m).  All  other  coherence  between  current 
meter  pairs  at  this  frequency  and  at  all  higher  frequencies  is  insignifi¬ 
cant. 

Table  3  shows  that  in  the  low-frequency  range  (which  includes  the  storm 
frequency)  high  coherence  is  found  between  some  current  meter  pairs  but  not 
others.  The  current  at  C  (NBA569)  is  well  correlated  with  the  current  at 
E  (NBA570),  for  example,  but  not  with  the  current  at  B  (NBA566).  The 
current  at  E  (NBA570)  is  wel 1 -correl ated  with  the  currents  at  both 
D  (NBA566)  and  F  (NBA567),  but  that  at  D  is  not  well  correlated  with  that 
at  F.  This  could  be  explained  by  E  being  up-stream  from  D  and  F  but  the 
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currents  at  D  and  F  belonging  to  different  eddies.  None  of  the  current 
meters  below  the  thermocline  is  wel  1 -correl ated  with  those  above,  except 
perhaps  NBA570  (E,  68m)  with  VACM355  (C,  14m). 


TABLE  3 

CURRENT-CURRENT  AND  CURRENT-WIND  COHERENCE  IN  LOW-FREQUENCY  RANGE 

(19  h  -  17  day  period)  the  asterisk  refers  to  range  (13-19  h) 


(B)  NBA565 

.72 

.23 

.52 

.36 

.63 

.12 

.21 

.29 

.67 

(D)  NBA566 

.12 

.42 

.47 

.73 

.29 

.27 

.23 

.45 

(F)  NBA567 

.36 

.56 

.78 

.24 

.15 

.21 

.46 

(A)  NBA568 

.69 

.72 

.13 

.15 

.29 

.80 

(C)  NBA569 

.75 

.23 

.11 

.33 

.36 

(E)  NBA570 

.45 

.26 

.34 

.46 

(C)  VACM355 

.44 

.20 

.95* 

.15 

(C)  VACM356 

.59 

.16* 

.15 

(C)VACf1357 

.10* 

.20 

(C)  NB3 

. 

0 

120 


SACLANTCEN  SR-75 


SUMMARY  AND  CONCLUSIONS 

The  data  collected  from  the  MILEX-80  experiment  have  been  presented  and 
i nterpreted. 

A_ description  of  the  experiment  is  provided  in  Ch.  2.  In  short,  six  ther¬ 
mistor  chains,  ten  current  meters,  and  a  meteo  chain  were  moored  by  means 
of  surface  floats  and  taut  nylon  rope  in  the  north  Tyrrhenian  Sea.  The 
instruments  were  moored  in  the  surface  mixed  layer  and  thermocline  and 
arranged  in  a  horizontal  distribution  to  cover  the  region  of  a  large  quasi¬ 
permanent  cold  patch  of  water  fanning  out  eastwards  from  the  Strait  of 
Bonifacio. 

The  cold  patch  has  been  interpreted  as  an  upwelling  caused  by  wind-stress 
curl  —  the  west  wind,  which  predominates  throughout  the  year,  being  con¬ 
sistently  stronger  at  the  eastern  end  of  the  Bonifacio  Strait  than  it  is 
off  the  east  of  Corsica  or  Sardinia  on  either  side  of  the  Strait.  This  up- 
welling,  commonly  referred  to  as  Ekman  pumping,  formed  the  central  theme  of 
the  report  and  set  the  context  for  the  various  spatial  and  temporal 
descriptions  of  the  data  occurring  throughout. 

The  existence  of  two  counterrotating  eddies,  one  cyclonic  underlying  the 
cold  patch  and  the  other  anti  cyclonic  in  the  warmer  water  to  the  south,  was 
confirmed  both  by  the  dynamic  topography  estimated  from  CTD  casts  and  from 
the_  measured  currents.  It  was  found  that  the  density  contribution  (i.e. 
estimated  dynamic  topogaphy)  to  the  eddies  accounted  for  approximately  only 
half  the  measured  speeds,  although  the  agreement  in  direction  was  quite 
good.  The  density  contribution  is  interpreted  as  a  baroclinic  adjustment 
to  forcing  by  the  wind-stress  curl.  The  remainder  of  the  flow  is  con¬ 
sidered  to  be  barotropic  but  not  necessarily  caused  by  local  forcing.  It 
might,  for  example,  be  an  instability  on  the  mean  circulation  in  the 
Tyrrhenian  Sea  as  it  flows  northward  up  the  west  coast  of  Italy,  turns  west 
and  south  under  the  influence  of  the  northern  barrier  to  the  basin,  and 
then  turns  southwards  down  the  west  of  Corsica  and  Sardinia. 

A  model  after  De  Szoeke  <34>  was  suggested  for  the  mixed  layer,  in  which 
the  mixed  layer  represents  a  balance  between  the  upward  motion  of  the  layer 
base  due  to  Ekman  pumping  and  the  downward  progression  of  the  interface  due 
to  mixing.  Then,  taking  a  wind-stress  curl  equivalent  to  4  m/s  over  100  km 
(a  value  that  appears  consistent  with  historical  data  presented  in 
May  <10>),  and  taking  an  entrainment  rate  (after  Kraus  &  Turner  <6>)  due  to 
a  wind  speed  of  8  m/s,  it  was  shown  that  a  steady  state  would  prevail  after 
7  days.  In  this  case  the  mixed-layer  depth  would  be  15  m  and  the  Ekman 
pumping  velocity  would  be  2.1  m/day.  The  time  scale  is  similar  to  the 
storm  period  and  the  depth  scale  close  to  the  measured  depth  of  15  to  20  m 
for  the  mixed  layer  over  the  cold  patch. 

Surface  heat  fluxes  were  estimated  for  both  the  cold  water  patch  and  the 
warmer  water  to  the  south  using  the  Bowen  ratio  technique.  This  involved 
taking  an  empirical  form  for  the  Bowen  ratio  dependent  on  air/sea  tem¬ 
perature  difference  and  estimating  from  this  both  the  evaporative  heat  loss 
and  the  water  vapour  correction  term  in  the  back-radiation  formula.  Cloud- 
cover  correction  factors  for  the  three-week  period  were  not  known  but 
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cloud-cover  statistics  for  the  two  areas  were  assumed  not  to  differ  signi¬ 
ficantly.  A  heat  flux  difference  for  the  two  areas  (north  &  south)  was 
thus  estimated  and,  by  comparison  with  the  north-south  difference  for  the 
heat  content  change  over  a  13-day  period,  a  vertical  velocity  difference 

-  w*5  between  north  and  south  was  estimated  and  found  to  be  5.9  m/day. 
By  considering  maximum  and  minimum  cloud-cover  conditions  this  implies  an 
upwelling  to  the  north  in  the  range  from  2.1  to  3.8  m/day  and  downwelling 
to  the  south  in  the  range  from  3.9  to  2.2  m/day.  The  wind-stress  curl 
necessary  to  produce  such  up-  and  down-welling  is  consistent  with  the 
historical  data. 

Another  possible  technique  for  estimating  wind-stress  curl  emerged  from 
Ch.  5,  on  spectral  analysis.  There  it  was  argued  that  the  peak  in  tem¬ 
perature  variance  spectra  occurring  at  the  inertial  frequency  could  be 
attributed  entirely  to  inertial  pulsations  in  the  divergence  of  the  Ekman 
mass  flux.  Assuming  a  simple  model  for  the  inertial  oscillations  arising 
from  the  sudden  application  of  a  constant  wind  stress,  wind-stress  curl  was 
estimated  as  a  function  of  temperature  variance,  the  inertial  frequency, 
and  the  mean  maximum  thermocline  gradient.  Although  involving  a  number  of 
assumptions  and  approximations,  the  technique,  rather  surprisingly,  gave 
values  of  the  wind-stress  curl  and  upwelling  in  good  agreement  both  with 
estimates  from  the  heat  balance  and  with  the  historical  data.  The  mean  up¬ 
welling  over  the  three-week  period  for  position  B,  for  example,  was  esti¬ 
mated  to  be  2.25  m/day. 

Similar  arguments  invoked  to  explain. the  inertial  peak  in  the  anticlockwise 
rotary  spectra  for  the  current-meter  data  were  less  successful. 
Nevertheless,  it  was  demonstrated  that  this  energy  could  arise  through  a 
resonant  forcing  of  the  subsurface  (below-thermocline)  layer  by  inertial 
pulsations  in  the  hydrostatic  pressure  gradient.  These  in  turn,  it  was 
argued,  are  a  consequence  of  the  pulsations  in  Ekman  mass  flux  divergence. 

The  following  more  general  statements  can  be  made  as  a  result  of  the 
MILEX-80  experiment: 

a)  The  spatial  variability  of  the  wind  field  over  the  Mediterranean 
strongly  influences  the  depth  of  the  mixed  layer  and  the  geostrophic 
current  field. 

b)  In  such  conditions  of  high  spatial  variability  in  the  wind  field 
one-dimensional  mixed-layer  models  cannot  work  well. 

c)  One-dimensional  mixed-layer  models  could  be  generalized  to  three 
dimensions  and  used  if  the  wind  field  were  known. 

d)  Areas  of  high  wind  variability  can  be  identified  both  by  use  of 
satellite  infrared  images  and  by  analysis  of  historical  data. 

e)  Many  such  areas  exist  in  the  Mediterranean. 

f)  Knowing  the  wind  field  would  allow  the  use  of  mixed-layer  models  in 
predicting  the  evolution  of  surface-layer  variability.  In  order  to  pro¬ 
perly  predict  surface-layer  variability,  however,  an  ocean  circulation 
model  is  required,  since  only  in  this  way  can  the  important  effects  of 
horizontal  advection  be  accounted  for. 
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The  area  of  MILEX-80  has  turned  out  to  be  too  complex  for  the  proper 
testing  or  application  of  one-dimensional  mixed-layer  models.  However, 
data  collected  during  a  second  mixed-layer  experiment,  MILEX-82,  performed 
in  the  open  Western  Mediterranean,  half  way  between  the  Balearic  Islands 
and  Sardinia,  should  prove  more  suitable  for  one-dimensional  model  tuning. 
This  is  due  to  the  distance  of  the  area  from  coastal  boundaries,  the  deep 
water  (over  2500  m),  the  flat  bottom,  and  the  low  spatial  variability  in 
the  wind-stress  field.  Furthermore,  satellite  infrared  images  studied 
during  the  MILEX-82  experiment  did  not  reveal  any  strong  gradients  in  the 
surface  temperature,  as  can  be  found  permanently  in  the  MILEX-80  zone. 

Further  research  will  involve  using  MILEX-82  data  to  tune  one-dimensional 
models  to  establish  appropriate  values  of  the  empirical  constants,  genera¬ 
lizing  these  models  to  three-dimensions,  and  finally  modelling  the  surface 
layer  in  the  MILEX-80  zone  by  assuming  the  spatial  variations  in  the  wind 
field  given  by  the  available  historical  data. 
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APPENDIX  A 


A  VERTICAL  VORTICITY  BALANCE  IN  NON-DIMENSIONAL  FORM 


Neglecting  molecular  viscosity,  a  vorticity  balance  can  be  written  in  the 
form  (<A.1>  p.  38) 


where  ^  is  the  absolute  vorticity  (i.e.  vorticity  of  fluid  relative  to 
the  earth  plus  vorticity  of  fluid  due  to  earth's  rotation),  ^  the  fluid 
velocity,  p  the  density,  and  A|-|  and  Ay  are  eddy  viscosity  coef¬ 
ficients  in  the  horizontal  and  vertical  respectively.  The  usual 
righthanded  cartesian  coordinate  system  (x,y,z)  is  assumed  with  z  directed 
(+ve)  upwards  normal  to  the  earth's  surface. 

The  second  term  of  Eq.  A1  represents  vortex  stretching  and  tilting,  the 
third  term  baroclinic  production  of  vorticity,  and  the  fourth  term  dissipa¬ 
tion  due  to  the  turbulent  Reynolds  stresses. 

Since  we  are  primarily  interested  in  the  effects  of  horizontal  stresses 
applied  to  the  sea  surface  by  the  wind  we  consider  only  the  vertical  com¬ 
ponent  of  Eq.  A1  given  by 


(Eq.  A2) 


with  w  =  (5,  n,  ?)  being  the  relative  fluid  vorticity  and  f  a  vertical 
vector  for  the  Coriolis  parameter  f  .  Dimensional  variables  Tn  Eq.  A2  are 
changed  to  non-dimensional  variables  by  the  transformations 


(x,  y)  [L]  (x,  y) 

z  >  [D]  z 

t  ^  [L]  t 


U 

(u,  v)  ^  [U]  (u,  v) 

w  [M]  w 


(Eq.  A3) 


L 
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P  [po]  P 

p  ->•  [po  ^9]  P^  [pq  ^  U  L]  p  , 

where  all  terms  in  square  brackets  are  dimensional,  with  L  a  horizontal 
length  scale,  D  a  vertical  length  scale,  U  a  horizontal  velocity  scale, 
and  Pq  a  density  scale. 

Substituting  the  transformations  of  Eq.  A3  into  Eq.  A2  and  multiplying 
throughout  by  L/fU  leads  to  the  non-dimensional  equation 


Rq  +  3v  +  Ro(^*V)w  +  -  2 _  r 

Dt  “  3z  p2  L 

-  Eh(VH  V§  u^)  -  Ev  ^Vh  X 


3p  3p 
3x  3y 


iP 

9y 


9£' 

3X  . 


3z2 


=  0 


(Eq.  A4) 


where  Rq  ,  Eh  and  Ey  are,  respectively,  the  non-dimensional  Rossby 
number,  horizontal  Ekman  number,  and  vertical  Ekman  number.  These  are 
defined  by 


_  ^  (Eq-  A5) 

"  fL2 


The  non-dimensional  parameter  3  in  Eq.  A4  is  expressed  by 


3 


(Eq.  A6) 


where  3o  is  the  south-north  gradient  in  the  linear  approximation  to  the 
variation  of  f  with  latitude  f  =  ^q  +  3oy  . 
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